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in recent years a considerable amount of theoretical, experimental, and computational work in 
the development and application of techniques for accurate measurement of microwave antennas has 
been successfully completed at the National Bureau of Standards (and work is continuing). This paper 
presents and extends the basic plane-wave scattering-matrix formalism and presents new generalized 
or adjoint reciprocity relations for antennas. The P WSM formalism is eminently suitable lor the formu- 
lation and solution of problems involving interactions at arbitrary distances and for the expression of 
conventional asymptotic quantities, such as gain, effective area, and polarization. It has in particular 
enabled derivation of two new techniques that permit accurate, "probe-corrected" antenna measure- 
ments at greatly reduced distances: (1) by deconvolution of transverse scanning data, taken with d <^ da 
(where cfo=a 2 /2A) and (2) by extrapolation of received signal observed as a function of distance <L 
with d ~ (Ir. These techniques basically determine the scalar product, C, of two vectors characteristic 
respectively of the transmitting and the receiving antennas. Formulas for utilization of (7-data, taking 
full account of polarization characteristics and not requiring reciprocal antennas, are given for (a) 
one-unknown-antenna, (b) generalized two-identical-antenna, and (c) generalized three-antenna measure- 
ment techniques. 
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Some Conventions and Notation Used Consistently 

(a) Conventions 
1. Complex numbers are sometimes called phasors; they are not called vectors. 



2. Bold-face symbols denote vectors or dyadics denned in "ordinary" space or in wavenumber 
space. Components may be complex numbers. 

3. Scalar and vector products of two vectors are denoted by A • B and AxB, respectively. 
The scalar product of three vectors taken in the cyclic order A, B, C is denoted by [ABC]. 

4. A superposed bar denotes the complex conjugate. 

5. The magnitude of a complex number z is denoted by \z\. 

6. The squared magnitude of a vector V is defined by V • V and denoted by | Y | 2 . 

7. The "square" of a vector V is defined by V • V and denoted by V 2 . Example: k 2 = o) 2 fie (see 
the list of symbols following). 

8. "Transverse" means perpendicular to the z axis unless otherwise indicated. 

9. "On-axis" refers to the z axis of coordinates, not to an axis possibly suggested by antenna 

geometry. 
10. The (suppressed) time dependence is exp (— io)t). 
(b) Roman Letters 
a Q : Complete vectorial spectrum for E of incident field ( p. 14) . 
A Q : Transverse part of a q (p. 15). 

ao: Incident wave-amplitude in antenna feed transmission line or waveguide (p. 9) . 
h Q : Complete vectorial spectrum for E of scattered or radiated field ( p. 14) . 
B^: Transverse part of h Q ( p. 15) . 

bo: Emergent wave-amplitude in antenna feed transmission line or waveguide. 
c?K: Symbolizes surface element in double integrals in k x , k y space. 
dK: Symbolizes surface element in double integrals in x, y space. 
E: "Electric field" (complex representation). 
Efi Transverse part of E. 
e x , e y , e z : Fixed, orthogonal, right-handed system of unit vectors. 
e||, ej_, efei Orthogonal, right-handed system of unit vectors tied to k(p. 13). 

G Q (K): Power-gain function evaluated in the direction of k; q= 1 or 2 implies k = k + or k~, 
respectively (p. 23). 
H: "Magnetic field" (complex representation). 
H*: T ransv erse part of H. 
i: V^ 

k: Propagation vector; components k x , k y , k z . 
k ± : Propagation vector with z-component equal to ±7. 
K: Transverse part of k (K is chosen real in this work). 
A : Vk • k= qj VjJie (a real quantity in this work). 
K: VK^K 
m, n: Index taking on values 1 , 2 and indicating association with unit vectors k x , K 2 (p. ). 

p, q: Index taking on values 1, 2 and indicating association with regions to the "right" and 
to the "left" of an antenna or scatterer (p. 12). 
r: Position vector (a real vector); components x, y, z. 
R : Transverse part of r. 
r: Magnitude of r. 
R : Magnitude of R. 
Sog : Transverse vectorial receiving characteristic (p. 17) . 
so<? : Complementary receiving characteristic ( p. 25) . 
Sgo : Transverse vectorial transmitting characteristic ( p. 16) . 
s 9 o : Complete transmitting characteristic (p. 23) . 
Wq(K) : Polarization index for incident plane-waves (p. 25). 

Yq : Ve//x, wave admittance for simple plane-waves in medium with parameters 6, (jl. 
(c) Greek Letters 

y: V/.' 2 — K 2 , taken positive when K < k :, positive imaginary when K > k; k z — ±7 (p. ). 
d(kjc) : Dirac delta "function." 
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8(K) : Abbreviation for 8 ( kj-) o (k y ) . 

€*, to: Permittivity of homogeneous, isotropic, dissipationless medium. 

77,,: Characteristic admittance for the propagated mode in waveguide feed (p. 9). 

tj , : z-component of wave-admittance for TM plane-waves in space; 171 = coe/y = Yok/y, 

and for ^ 0*Z 7r, = y /|cos 0| (p. 12). 

iqr. z-component of wave-admittance for TE plane-waves in space; r)2 = yl (co/jl) = Yoy/k, 
and forO^ 0*5: tt, = Fo|cos 0j(p. 12). 

0: Polar angle in spherical polar coordinates (p. 13). 
K\: Unit vector = K//\(p. 11). 
k> : Unit vector = e z X ici (p. 1 1) . 

H , fio : Permeability of homogeneous, isotropic, dissipationless medium. 

p<> 9 : Polarization index for receiving characteristics ( p. 25) . 

Pqo : Polarization index for transmitting characteristics ( p. 24) . 

o-g(K) : Effective area for reception (p. 24) ; g= 1 or 2 implies k=k or k = k f , respectively. 

4> : Azimuthal angle of plane or spherical polar coordinates ( p. 13) . 

<*) : Angular velocity as in the suppressed time factor exp (— iwt). 

Introduction 

In section 1 we give, quite thoroughly as far as definitions and notation are concerned, a formu- 
lation of the plane-wave scattering-matrix for antennas and scatterers. ' In previous publications 
[1, 2], 2 only a "one-side" matrix description of antennas, using one spatial reference plane, was 
given explicitly. To provide a more complete exposition, we give here the formulation of the 
complete scattering matrix obtained by enclosing the antenna between two planes and considering 
incident and emergent plane waves on both planes. 

The plane-wave scattering matrix is put forth as a good basis for advanced antenna theory, 
especially such theory related to antenna measurement techniques. While the results given may be 
persuasive, numerous analytical examples and applications which illustrate and extend the theory 
are available — and are not included. These are intended for a more comprehensive publication. 
Apart from appendices A and B, relatively little theory is developed in this paper. However, logical 
interrelations are indicated, and more detail is given when it is a question of relating new concepts 
to more familiar concepts and to practicable antenna measurement techniques. 

An "acoustics translation" of much of the present work has been published [3J. Because of 
the relative simplicity of the acoustic wave-fields involved, this may be found helpful in illuminating 
the electromagnetic case. 

The scattering matrix formulation properly includes basic expressions for power transfer and 
for reciprocity. 

The reciprocity relations are stated in a generalized form, using the concept of mutually ad- 
joint antennas. The concept of generalized reciprocity is not in itself new [4, 5], but apparently it 
has not previously been formulated for antennas and scatterers (see, however, remark following eq. 
(1.6-21)). The generalized or adjoint reciprocity relations have found substantial application in 
research establishing the foundations of the extrapolation technique [6, 7]. They are used in this 
paper in the formulation of possible new antenna measurement techniques, predicated upon the 
physical realization of mutually adjoint antennas. 

In section 2 the plane-wave scattering matrix approach is used to obtain a complete and general 
solution to the problem of coupled antennas. In spite of their formal appearance, the general re- 
sults obtained represent the heart of the present theory, and provide a fruitful and reliable basis 
for additional results (including those reported here). Two of these results are the deconvolution and 
extrapolation techniques, labeled (1) and (2) in the Abstract and outlined in section 3. These techni- 



1 Since an antenna is in general a scattering object, and a passive antenna externally is merely a scattering object, the consideration of scattering is included in 
the full consideration of antennas and need not always be mentioned explicitly. 

2 Figures in brackets indicate the literature references at the end of this paper. 



ques basically determine values of the scalar product, called the coupling product, of two two-com- 
ponent vectors characteristic of the two antennas involved. (This statement also applies to the conve- 
tional far-field antenna measurement methods, provided due attention is paid to polarization charac- 
teristics.) 

Utilization of coupling-product data is discussed as a separate topic. Here the analytical prob- 
lem is primarily geometric and algebraic, with the exact form depending upon what is considered 
known a priori and what information is sought. Three classes of antenna measurement situations 
[labeled (a), (b), and (c) in the Abstract] are discussed in section 4. The order of listing and dis- 
cussion is roughly that of decreasing a priori information and increasing complexity. 

An increasing body of experimental results involving various combinations of the several 
techniques identified above may be found in the literature [8-11, 54]. Moreover, the error analysis 
required for determination of accuracy in concrete measurement situations is approaching com- 
pletion (Kanda [ 12] , Yaghjian [ 13] , Newell [14]). 

The combination of the techniques labeled (1) and (a) amounts to a technique for correction of 
near-field antenna measurements made with an arbitrary but known measuring antenna. The ability 
to obtain the true radiated spectrum of an unknown antenna, fully corrected for the effects of the 
measuring antenna, incidentally implies the ability to obtain corresponding true values of E and H 
in the near field, similarly fully corrected. Frequently, and in particular in the following paragraphs, 
the measuring antenna will be referred to as a "probe." 

The general subject of determination of far-field antenna patterns from near field data is 
surveyed and an extensive bibliography is given in a recent paper by Johnson et al. [15]. This paper 
should be consulted for an overview of the subject. Here we mention specifically only certain earlier 
work in which the "probe-correction problem" was considered or which represented steps leading 
to the eventual simple, rigorous, and general deconvolution solution of the problem. 

Woonton. in 1953 [16], obtained an integral expression for the near- field response of a linear 
(=thin wire) antenna and discussed probe effects qualitatively. Woonton stated that the problem 
had not been critically discussed previously. Dayhoff (1956) [17], using scalar waves, plane-wave 
spectrum analysis, and reciprocity, introduced a version of the very important transmission integral. 
(Dayhoff used the transmission integral to obtain an approximate solution of the diffraction cor- 
rection problem in microwave interferometry. A rigorous and more general version of this solution 
was presented by Kerns in 1957 [8|.) Brown [19] (1958), using plane-wave spectrum analysis and 
reciprocity [20], obtained a version of the transmission integral and used it to give an approximate 
analysis of probe effects, assuming simple, known data for both antennas involved. More work along 
this line (limited to two-dimensions) was done by Jull [21, 22]. In 1961 Brown and Jull [23] gave a 
rigorous and general solution to the probe correction problem for the two-dimensional case using 
cylindrical wave functions. The use of two-dimensional solutions in certain three-dimensional 
problems has been suggested, and was studied experimentally by Martin [24], but is not valid for 
any three-dimensional problems. The proper extension to three dimensions in spherical or cylin- 
drical coordinates is far from trivial. See Jensen [25], Leach and Paris [26], and Wacker [27]. 

The key to the present solution to the probe-correction problem is the use of plane-wave 
analysis and rectangular coordinates and the recognition that planar scanning would permit 
solution of the integral equation presented by the transmission integral by Fourier inversion 15 
(better called deconvolution in the existing context). This solution was presented in 1963 by Kerns 
[28], and again, including application of a two-dimensional, spatial sampling theorem, in 1967, 
at a University of Colorado Advanced Electromagnetic Theory Summer Course. Archival publica- 
tion, accompanied by substantial experimental application and verification by Baird et al. [8] , 
was accomplished in 1970 [2]. The transmission integral used, though similar to that derived by 
Brown, was actually obtained by methods used in Kerns and Dayhoff [1], where it was derived 
without recourse to reciprocity and was explicitly recognized as the first term in an infinite series 
of interaction terms. (See the remarks at the end of appendix A, below.) 



These matters are adequately discussed later in this paper. 



The Kerns and Dayhoff paper, not originally considered to be in the domain of antenna theory 
by its authors, has served as an important base for much of the material embodied in the present 
and related papers and in other work as yet unpublished. 
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Plane-Wave Scattering-Matrix for Antennas and Scatterers 
1.1. Representation of Fields on S , Definition of a and b 



Let us eonsider the antenna system shown schematically in figure 1. We choose a (mathe- 
matical) terminal surface So in the waveguide feed and define a supplementary surface S«, such 
that Sa + So forms a closed surface enclosing the source or detector associated with the antenna. 
The surface S a coincides with shielding, which is required to make the problem well defined 
(both experimentally and theoretically!). As an important measure of simplicity, we consider only 
the case of a single waveguide feed, supporting just one propagated mode. We employ conventional 
phasor wave amplitudes do and b {) for the incident and emergent traveling wave components at 
So. These wave amplitudes are fully defined by the following four equations. 

The tangential (= transverse) components of E and H on So are given by 



E f = (a + 6o)eo(r), 

Ho/=7}o(«o — Mh<>(r), i (r on So) 



(1.1-la) 
(1.1-lb) 



where eo(r) and h<>(r) are real basis fields for the mode involved, subject to the impedance 
normalization 



ho(r) = 7),, l T7u,iio X e (r) 



(1.1-2) 







Oi 



(z = z 2 <o) 



Zj>0) 



FIGURE 1. Some notation for plane-wave scattering-matrix description of antennas. 



Arrows indicate the association of as and 6's with incident and emergent waves respectively. Antenna-system representation is symbolic: no particular size, 
shape, symmetry, type, orientation, or position is implied. In any concrete case, position and orientation of antenna relative to coordinate system must be established. 
as with the aid of fiducial marks on the antenna. 



and to the power normalization 

[eoh„no]dS=l. (1.1-3) 



/ 



So 



Here the integrand is the scalar triple product; no is the unit normal vector, drawn inward with 
respect to the antenna; 170 is the characteristic admittance, and 77 w is the wave admittance for the 
mode involved. Equation (3) establishes peak- value normalization for a ( > and bo, so that net time- 
average power input to the antenna at So is given by 

If— 1 

Po = ^Re EofXHo, -n dS=^T7o( | a | 2 - | 60 | 2 ), (1.1-4) 

z J s * 

where Re denotes that the real part is to be taken, the superposed bar denotes the complex con- 
jugate, and the vertical bars denote absolute values. Remarks: (a) The impedance normalization 
shown above is more flexible than that used previously [1, 2]. It allows one to choose whatever 
characteristic impedance or admittance is deemed familiar or convenient for purposes of measure- 
ment-related calculations. Examples are, for waveguide, 170=1 or r)o=r) w \ and, for coaxial line, 
the conventional characteristic admittance ( = 2tt V€//x,/ln(6/a) in conventional notation), (b) 
Power normalization differs slightly from that used previously: Powers of the factor (2tt) have 
been redistributed, (c) For a detailed discussion of the material in this subsection and for the small 
but essential amount of microwave network theory needed in measurement-related calculations, 
see [29] or the especially prepared report [30]. For more elementary discussion see also [31]. 

1.2. Representation of Fields in Space; Definition of a q (m f K) and b u (m, K) 

We choose a rectangular coordinate system Oxyz (with unit vectors e^, e y , and e z ) so that the 
considered antenna system may be confined entirely to the space between the (mathematical) sur- 
faces F\ and F> at z = Z\ ^0 and at z = z-> ^ 0, as suggested in figure 1. The electromagnetic fields in 
the regions to the "right" and to the "left" are to be represented as superpositions of plane-wave 
solutions of Maxwell's equations. This type of representation is well known (see e.g., [32)), at least 
for solutions of the scalar Helmholtz equation; an appropriate generalization to the (vector) elec- 
tromagnetic field, though often shunned, offers no particular difficulty. 

The electromagnetic field in the regions under consideration satisfies Maxwell's equations in 
the form 

VxE-iw^tH, VxH = -z"coeE, (1.2-1) 

where /it, e are constant real scalars representing respectively the permeability and the permit- 
tivity of the medium, and exp (— ixot) time dependence is assumed. We derive our basis fields from 
the general plane wave 

E = Texp(*k-r), \ 

> (1.2-2a) 

H=(co/a)- 1 kxTexp (ik-r)J 

which is a solution of (1) for any propagation vector k such that & 2 = k-k=o) 2 /A€ and any vector 
T (independent of position r) satisfying the transversality relation 

k-T=0. (1.2-2b) 

In spite of this occurrence of "transversality," in what follows the term "transverse" will mean 
transverse with respect to the z direction unless otherwise specified. 

The propagation vector will be regarded as a function of its transverse components k x , k y 
(which are chosen real); the z-component is thus 

k z = ±y, (1.2-3a) 
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where y 2 = A; 2 — k%— k' y . It will be convenient to denote the transverse part of the propagation 
vector by K, so that K= k x &x + k y e y and 



y- 



(tf-K 2 ) 1 ' 2 . 



(1.2-3b) 



Since k x , ky must be allowed to vary independently in the range (— o°, o°) 5 real and imaginary values 
of y will occur, y will be taken positive for K~ < /r, positive imaginary for K 1 > k 2 . Superscripts 
4 N-" and "— " will be used when it is desired to indicate the choice of sign associated with k z . When 
y is real, the exponentials exp (ik + -r) and exp (ik -r) respectively represent simple plane waves 
traveling into the +z and — z hemispheres. When y is imaginary, the exponentials represent in- 
homogeneous plane waves with propagation of phase in the transverse directions and exponential 
attenuation of amplitude ("evanescence") in the + and — z directions, respectively. 

In virtue of the relation k*T=0, (2a) yields just two linearly independent fields, hence just 
two basis fields, for any given k. The appropriate polarizations for the basis fields are those with 
the electric vectors parallel or perpendicular to the plane of k and e 2 , which is the plane of in- 
cidence for a wave incident on any plane z = const. This choice of polarizations yields "trans- 
verse magnetic" and "transverse electric" waves; the same choice of polarizations, usually labelled 
"I?//" and "2?j_", simplifies the derivation of FresneFs equations in optics or electromagnetic 
theory. 

In order to set up the basis fields in the desired form, we require the transverse unit vectors 



K\ 



K//C, K 2 = e z X K u 



(1.2-4) 



which are respectively in and perpendicular to the plane of k and e z . This part of the notation is 
illustrated in figure 2; K\ and k> may be identified as radial and tangential unit vectors, as often 




FIGURE 2. Illustrating k, K, k. and #o. 



associated with polar coordinates in a plane. As a temporary abbreviation 4 we put // ! =exp (ik ± -r)/ 
(27r). For the E ]{ (or TM) components we put T= K\ TKy~ l e z and obtain from (2a) 



E, = [ki -hKy l e z \u ± , 
H 1 ± = ±T7ie 2 X Ki^ ± , 



(1.2-5) 



4 The factor ll(2n) in the definition of u represents a change in normalization consistent with that noted for a and bo in the preceding subsection. 



11 



where 171 = coe/y. For the E ± (or TE) components we take T = k> and obtain from (2a) 

\ (1-2-6) 

H^t+^X^ + ^l-'e,]^,) 

where j]> — y\ (w/x). Among other similarities it may be observed that 171 , 172 are wave-admittances 
that correspond closely to the wave-admittances encountered in the theory of rectangular wave- 
guide. Equations (5) and (6) furnish the desired basis fields; somewhat arbitrarily, we have chosen 
to make the expressions for transverse E as simple as possible. The normalization and orthogonality 
properties of the basis fields are of course implicit in the expressions themselves. (This is an inter- 
esting contrast to (1.1-1), where the field patterns are implicit and the normalizations explicit.) 

Let us now examine the plane-wave representations of electromagnetic fields in the regions 
z ^ z\ > and z ^ z 2 < 0. We write 

E «( r )= f° [ 2 EM m ' K)E^(K,r) + a g (m, K)Ej&(K,r)]d*,dfc y , (1.2-7a) 

J - x J m = 1 

H g (r)=[°° [ J) [M^K)H*(K,r) + a ff (m,K)H^K,r)]rffr^ y . (1.2-7b) 

J -00 J m=l 

Here the index q takes on the values 1 and 2 and identifies quantities respectively associated with 
the regions to the "right" and to the "left" of the system considered; the upper and the lower 
superscript signs are associated with q= 1 and q=2, respectively; and b q (m, K) and a q {m, K) 
are scalar spectral-density functions for outgoing and incoming waves, respectively. The electro- 
magnetic fields given by (7) will satisfy Maxwell's equations provided that the necessary differentia- 
tions can be taken under the integral signs. 

Now, as will be shown in a moment, a knowledge of the transverse components of E r/ and H Q 
(in a single plane, in fact) is sufficient to determine a Q (m, K) and b q {m, K); and hence, by (7), the 
entire electromagnetic field in each of the regions considered. The z-components of the fields in 
(7) are, strictly speaking, redundant. The transverse components are both necessary and suf- 
ficient for the expression of normal energy-flux and continuity conditions across a transverse plane. 
We find, in fact, that the inclusion of z-components is sometimes convenient and sometimes not. 
The following equations illustrate the latter case. 

For the transverse components of E and H in the regions z 5* z\ ^ and z ^ z 2 ^ 0, we find from 
(5), (6) and (7) the Fourier integral representations 

E#(r)=^[]£ [b q (m,K)eiy\*\+a q (m,K)e-iy\*q Km e iKB dK % (1.2-8a) 

H ff ,(r)=^ J J [b q (m, K) e *W a q (m % K)e~* w ]n q X KmJ\ m (K)e iK R dK, (1.2-8b) 

where 171= coe/y, r]> = y\ (co/jl) , and r = R-fze,. Here and in subsequent expressions of this 
type summation over the values 1 and 2 of the polarization index m and integration over the in- 
finite kjc, k y plane is to be understood. The role of the index q is as described under (7); the use of 
\z\ and — \z\ assures proper phase variation for outgoing and incoming waves, respectively; and 
n,/ is the outward normal unit vector on F q : ni = e^, n 2 = — e^(i.e., "outward" means with respect 
to the slab between the surfaces Fiand F 2 ). The functions a q (m, K) and b q (m, K) maybe regarded 
as modal terminal variables for the continuous spectrum; a mode is identified by a triplet of values 
(m, k x , k y ) and a direction (right ward or leftward). 
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Explicit expressions for the spectral density functions b lf {m, K) and a q (m, K) may be found 
from the Fourier inversion <>f the above equations: 

b q (m, W) = 6 -^-k„, • [ [E,(R,z)+7,^H g (R,2) X n Q ] e -' KR dR, (1.2-9a) 

a Q (m, K) = e -^— Km ■ f[E,(R,z)-T);;H,(R,z)Xn,]e"' K '>. (1.2-9b) 

477 J 

Here the integrations are to be taken over the entire x, y plane — as required by the Fourier in- 
version—for any suitable fixed value of z (for q = 1, z 5= Z\\ for q = 2, z ^ z>i). 

The spectral functions b q (m, K) and a q (m, K) are independent of z, although this may not 
be immediately apparent in (9), and the values of the functions are referred to the plane z = 0. If 
desired, the phases and amplitudes could be referred to other reference planes S (j at z= Z q , say. 
However, the choice implicitly made, Z\ = Z> = 0, is convenient at least for present purposes. 

It is convenient at this point to introduce a number of definitions for future reference when 
and as needed. 

We observe that (7a) can be written 

E v (r)=-^- I [b,(K)e'>l^ + a,(K)e-'>l^l]^ K ll r/K. (1.2-10) 

Z7T J 




FIGURE 3. Unit vectors associated with k and K (for real kj. See also table 1. 
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Table 1. Unit Vectors Associated With K and k 



K: (A-. r , k y ) = (Kcos4>, Ks'mcf)) 



Ki = K//C = (A\ r e. r + k u e y )IK = cos cf> e x + sin e y 

Ki = e z X Ki = (— A\ v e.r + k. r e,j)/ K = —sin $ e^ + cos (/) e,, 

[KiK 2 C 2 ] = 1, Km(-K) =~Km(K) 



k: (A>, Zr y , A:*) = (A sin # cos </>, A sin sin c£, A cos 6?) 



ek — k/A = (/Cki + k z e z )/k = sin 6^ #ci H- cos 6 e z 
e || = K 2 X ca- = (A: 2 Ki — Ke 2 )/k = cos 6 Ki — sin e 2 

e ± = K2 

[e,|e x ejfc] = 1, e M (— k) = e„ (k), k z = ±y 



Notes: 

(a) K\ and #c 2 are respectively "radial" and "tangential" unit vectors associated with the plane polar coordinates 
for K. These unit vectors are always real. 

(b) In the propagating regime, where e N as well as ej_ is real, en and e_j_ may also be identified as the customary 6 
and unit- vectors of the spherical polar coordinates for k. The "middle" forms in the table help show what happens when 
K > A: en and e* become complex, but remain unit vectors in the sense en • e|| = e/, ■ e/,- = 1. 

(c) The unit vectors #o , k>, and en are not defined by the equations in the table at the singular points K = and 
6 — 0, 77. For k= Ae 2 , one may, e.g., determine a consistent set by choosing K\ — e x . (cf. use of (4.1-4)). 

Here h q and a,, the "complete vectorial spectra" for the outgoing and incoming plane-wave com- 
ponents of E 9 , respectively, are given by 

MK) = 6,(1,K) (Ki+Ky- l e z ) + b q (2, K)k 2 , (1.2-lla) 

Bq(K) = a q (l, K) (kj ± Ky-*e z ) + a q (2, K) k 2 . (1.2-llb) 

(The association of the upper and lower signs with q = 1 and q = 2, respectively, is continued from 
(7).) Alternatively, b q and a q may be expressed as follows: 

MK) = 6,(1, K) [(-)*- 1 */?] e (l (k) + 6,(2, K)cl(K), (1.2-12a) 

a 9 (K) = a,(l, K) [ Rft/y] e„ (k) + a q (2, K) e ± (K), (1.2-12b) 

where e |j (k) is denned as k 2 X k/A:, and e ± (K) is a suggestive alternative notation for #c 2 . (The 
factors in brackets in (12) are introduced as a convenient, explicit way of helping to keep signs 
straight.) The complete vectorial spectra are of interest to us primarily in the propagating regime, 
where ey, as well as ej_, is real. 

All the unit vectors associated with k and K are pictured in figure 3 and fully identified and 
related in table 1. 

The complete vectorial spectra must and do satisfy the transversality relations 

k • b 9 =0, k- • a 9 =0, (1.2-13) 

which indeed are expressions of the transversality of the basis fields (5) and (6). 
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We next observe that (8a) can he written 

Etf(r)=^:J [IMK)e'>M + A 7 (K)e-<>Mk' K ",/K, (1.2-14) 

where B f/ and A f/ , the "transverse vectorial spectra'' for the outgoing and incoming plane-wave 
components of E<#, respectively, are given by 

B 9 (K) = ^ b q (m, K)k w , (1.2-15a) 

m 

Ag(K) = ^ a </(™> K)Km. (1.2-15b) 

m 

Further we note that, given (15), z-components can be recovered with the aid of the transversality 
relations (13) or simply by inspection of (11). 

Remarks: (a) One may observe a switch — from scalar spectra and vector waves in (7) and (8) 
to vector spectra and scalar waves in (10) and (14). Equation (7) remains fundamental, (b) The use 
of the lower case letters a and b to denote complete spectra and the use of the capitals A and B to 
denote the corresponding transverse parts should be noted. This parallels our use of the letters 
r, k and R, K. 

The representations (7) and (8) afford a similar and essential resolution of any electro- 
magnetic field into two major parts: that represented by incoming waves and that represented 
by outgoing waves. The role of the former part is identified variously by terms such as incident, 
exciting, primary, or incoming; the latter, by terms such as induced, secondary, radiated, reradiated, 
scattered, or outgoing. The scattering matrix, defined in the next section, will he seen to be a way 
of specifying the linear transformation from the first part of the field to the second. 

The asymptotic relation between the far -field values of a scattered or radiated field and its 
spectrum is of essential interest and importance. This relation reveals the result of the interference 
among the waves of the continuous spectrum at large distances in any chosen direction of observa- 
tion. Normally this interference results (remarkably) in the well-known exp (ikr)/r variation with 
distance. This normal result does not apply to the spectrum of a plane wave, and it might not apply 
to the radiation or reradiation from a source of infinite size in one or more dimensions. For an 
antenna of finite size radiating into 3-dimensional space, we do have the asymptotic relations 
[33, p. 750J 

E£(r) ~ -ik | cos | B Q (Rk/r) e ikr /r, (1.2-16a) 

E r( r ) ik | cos | b q (Rk/r)e ikr /r. (1.2-16b) 

Here the superscript "r" refers to the radiated or reradiated component of the field and 6 is the 
polar angle of r with respect to the z axis. The first of the two equations relates to (14) and the 
second to (10); they differ only in the presence or absence of the z-component. Both equations are 
valid whether or not incident waves are also present in the field. Note that we have written RAr/r 
for K as the argument of the spectral functions. This expresses the fact that the vectors r and k 
involved must he parallel. In fact if we introduce spherical coordinates for r such that x = r sin 6 
cos (f), y= r sin sin $, z= r cos 0, we see that h q and B f/ are expressed as functions of the angular 
coordinates of r. Furthermore we note that here only real directions of propagation come into 
consideration, so that y is real and y = k | cos | is valid. 

To conclude this subsection, we give a useful, nonphysical, special result, which offers a 
sharp contrast to the type of results given just above. Namely, if an electromagnetic field is every- 
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where the (simple or evanescent) plane wave with E= a exp (iko • r)/27r, then from (9) the corre- 
sponding spectral functions are found to be 

ai(w,K)=b 2 (m,K) = a B i8(K-Ko),] 

(1.2-17) 
bi(m 9 K) =a 2 (m,K). = 0, J 

where a m — Km-a and 8(K — Ko) is an abbreviation for the delta-function product 8(k x — k , 
8(ky — k oy ). The vector 

A = a 1 Ki J ta 2 K 2 (1.2-18) 

will be known as the transverse spectral vector associated with the considered plane wave. 

1.3. Definition of Scattering Matrices for Antennas and Scatterers 

Let us first consider that a passive material structure, exhibiting linear electromagnetic 
behavior, is present in the region Z\ < z< z>. The scattering equations are written 

b ll (m,K)=y j f yS ff (ia,K;n,L)op(n,L)dL, (q= 1,2) (1.3-1) 



=2 f 2 



where, in addition to the summation and integration conventions noted following (1.2-8), we 
have summation over the index p, giving the contributions from the waves incident both from the 
"right" and from the "left." The processes described by the functions S u and S22 will be called 
backscattering; those described by Svz and S21 will be called transcattering. Essentially the same 
definitions and an example of (1) may be found in [1]. 

If the scattering object is also an antenna, the scattering matrix must include the transmitting 
and receiving characteristics, and the scattering equations are written 

fe = S()o«o + 2 I ^ So P (nX)ap{n,L)dL, (1.3-2a) 

p J n 

6,(/n,K)=S g0 (m,K)a + 2 J J) S 9P (m,K;n,L)a p (»,L)c/L. (1.3-2b) 



p y ^ /) 



Here #=1,2 and we have made use of the quantities 60 and a () , defined in (1.1-1). The quantity 
Soo represents "backscattering" observed at So in the feed waveguide, and the functions So 9 (m,K) 
and S q o{m^K) respectively represent the receiving and the transmitting characteristics of the 
antenna. (The quantities bearing the subscripts q or p = 2 in (2) represent the desired generaliza- 
tion of the antenna scattering matrix originally defined by Kerns and Dayhoff [1].) 

The definition of the antenna scattering matrix is now literally complete. We do not wish 
to belabor the generality of the definition, but we do call attention to the absence of restrictive as- 
sumptions. At best, its full significance can be made apparent only gradually. 

It will be advantageous to have the basic scattering equations expressed in vector-dyadic 
form. In addition to the transverse vectors B Q and A q defined in (1.2-15), we introduce the vectorial 
transmitting characteristic 

S q0 (K)=^S q o(m,K) Km, (1.3-3a) 
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the vectorial receiving characteristic 

S„„(K) = 2 S 0g (m,K) Km, (1.3-3b) 

in 

and the dyadic scattering characteristic 

S 9P (K,L) = ]T S qp (m,K;n,L) K m k n , (1.3-3c) 

m , n 

where the K m and k n are the unit vectors associated respectively with K and L. The scattering 
equations become 

6o = Sooa +2 f So„(K)-A„ (K)rfK (1.3-4a) 

B,(K) = S 90 (K)a + 2 f S w (K f L)-A„(L).rfL. (1.3-4b) 

The equations are now invariant with respect to choice of coordinates in the transverse plane. 
This facilitates discussion of "real world" cases, in which it is neither expedient nor necessary 
always to observe or produce 4 pure #o or #C2,field components. 

One may obtain an analysis of the scattering equations by considering the simplest modes of 
excitation: by a wave represented by do alone, and by individual spatial plane waves, represented 
by delta-function spectra. This procedure leads essentially to re-expressions of the basic defini- 
tions contained in the scattering equations. The circumstances of the resulting definitions or 
expressions are simple enough to suggest several more or less direct and conventional methods 
of measurement for the scattering-matrix elements. Let us in particular consider the definitions of 
S, /() (K) and Sog(K) — the scattering-matrix quantities that will receive the most attention in our 
work. 

Let us consider an antenna operating in its transmitting mode; that is, an antenna excited only 
by an incident wave in its waveguide feed and radiating into empty space. This elementary pattern 
of excitation is represented by a<> # 0, a q (m,K.) = 0. Under these conditions, the spectra radi- 
ated (to the left and to the right), normalized to unit ao, characterize the transmitting properties 
of the antenna; indeed, from (4b) 

B,(K) 

S 90 (K)=— . (1.3-5) 

ao 

Incident spatial waves being absent, (1.2— 9a) and (1.2-9b) together imply that B g may be related 
either to E or to H; in terms of E we have 

S 7 o(K)-^^ f E qt (R,z)e- iKR dR, (1.3-6) 

where z > Z\ or z < 22 must be in force. This gives us a Fourier transform definition of S 9 o, and, 
to the extent that E 7 //ao is known or measurable, a means of calculating S 9 o- The asymptotic 
relation (1.2-16a), applied in the present circumstances, may be written 

E^(r) ~ - iyS q o(Rklr)a e ikr lr. (1.3-7) 

Clearly this relation may be regarded as a formula for determining Sqo(K) in terms of observed 
asymptotic K qt lao. It in fact represents the basis of so-called direct methods of measurement 
of transmitting characteristics. 
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The combination of (6) and (7) incidentally furnishes a rigorous, vectorial form of the Fourier- 
transform relation between far- and near-fields. 

Next, consider an antenna excited solely by a spatial plane wave, incident on side q, say, and 
having E = a exp (ik*r)/27r. The pattern of excitation is accordingly represented by ao — 
and a q < (ra,K') = a )n 8 q > q 8 (K' — K), where a m =k,„ *a and the spectrum of the incident wave 
isfaund just as in (1.2-17). The scattering equations (2a) or (4a) now yield 

6o = So,(l, K)a 1 -hSo 9 (2,K)a 2 , 

(1.3-8) 
= S 0<? (KM. 

In other words, So q (m,K) denotes the receiving sensitivity — or receptivity — of the antenna 
to the Km-component of polarization of a plane wave incident on the antenna with direction of 
incidence specified by q and K. The normalization is to unit a m (which means unit delta-function 
spectrum amplitude). 

Equation (8) immediately suggests basic equations for direct measurement of receiving 
characteristics. Viz., 

&; = So g (K) -A', 

(1.3-9) 
&"=So g !(K) -A", 

where b' Q ,b" {) are observed and A', A" are known and linearly independent. The equation can, of 
course, be solved for So?; an algebraically identical problem is encountered in subsection 4.1. 
The following compact notation for the scattering equations, to be secured by introducing 
function vectors and making more use of matrix-algebraic concepts, provides a perspective quite 
different from that of the discussion just preceding. We first define the column matrices 



(1.3-10) 



in which a q and b q may themselves be regarded as column matrices representing the functions 
a (/ (m,K) and b q (m,\L). That is to say, the elements of these column matrices are labeled or 
indexed according to the values of m and K, and have the values a q (m,K.) and 6 r/ (m,K), re- 
spectively. The transformation from the entire set of incident waves to the entire set of emergent 
waves is now written 

\ /SoO Soi So2 \ / a °\ 

5i J = I Sio S u S l2 1 [ a, ] (1.3-11) 

\b-i/ \S 2 o S21 S22/ \o,2/ 

or, equivalently, after performing the indicated matrix multiplication, 

60 = Sooflo H" Soidi -t~So2fl2, 

b\ = Sioao + S\\a x -h Si2«2, (1.3-12) 

62 — S2OC0 + ^21^1 -h S 2 2«2- 
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Here the four kinds of products involved are defined by comparison with (2) or (4). This compact 
notation, in which the subscripts identify the three input-output reference surfaces, makes the 
overall structure of the equations more apparent and is a practical necessity for the demonstra- 
tion in the next section. The rules of matrix algebra apply: the S^ (as well as a q and b q ) correspond 
to column matrices; the So q9 to row matrices; and the S pq , to square matrices. 

We should notice the form that the scattering equations take in the absence of any scattering 
object (or conceivably in the presence of a non-scattering object). The free passage of waves is 
expressed by 

61 (m,K) = a 2 (m,K), b 2 (m,K) = a t {m 9 K). (1.3-13) 

(No propagation factors of the form exp (±iyd) appear here because the phase reference surfaces 
were chosen coincident at z = 0). The pertinent submatrix in (11) is 



(1.3-14) 



Here 1 denotes the identity transformation with elements S /W ,|8 (k x — L x ) 8 (k y — /,,) and denotes 
the zero transformation. One may regard a scattering object as producing a perturbation of the 
properties of free-space described in (13) or (14). From this viewpoint, what we may call true or 
bona fide transcattering is described by the operators S12 — 1 and S21 — 1. 

1 .4. Power Expressions 

The requisite expression for one-mode power transfer in the antenna feed waveguide was 
given in (1.1—4). Here we shall give the corresponding expressions for power transfer across the 
surfaces F\ and F>±. 

The time-average energy flux across the surface F q in the outward direction (the direction 
of n q ) is given by 



P g = jRef [EHnJdR, (1.4-1) 



where the integrand is the scalar triple product and the superposed bar denotes the complex con- 
jugate. We wish to evaluate P q in terms of the spectral density functions a q (m 9 \H) and 6 r/ (m,K). 
After some analysis, one obtains 

p * = \ \ S (|6«l 2 H«fll 2 h»<*- f £Re(W?m)dK, (1.4-2) 

Z jK<k m Ja>A- m 

noting that for the lossless medium y and r\ m are both real in the propagating region (K<k) 
and both imaginary in the evanescent region (K> k). The equation shows that power may be 
transferred by the coupling of incoming and outgoing evanescent modes having the same m and 
K. (This the essential mechanism of power transfer in a waveguide-below-cutofT attenuator.) 
We expect and assume that ordinarily this interaction will be negligible, as in the case of coupling 
of ordinary waveguide junctions; the separation required is measured in wavelengths, not Ray- 
leigh distances. Thus the basic expression (2) may ordinarily be abbreviated to 

P <i = ll E i\bq(m 9 K)\* - \a q (m 9 K)\ 2 ]7i m (K) dK. (1.4-3) 
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This may be used, e.g., to deduce conservation relations. In the absence of incident waves, (3) 
expresses the equivalent of "pattern integration" over a hemisphere. 

1.5. Reciprocal Relations — A Summary 

In this section we briefly state recently derived [34] generalized or adjoint reciprocity relations 
for antennas and scatterers. (The derivation is reproduced in appendix A.) These relations very 
readily adapt to the expression of ordinary reciprocity as a special case. We shall also comment on 
the question of "realizability" of mutually adjoint systems. First we must define mutually adjoint 
media and systems. 

We describe the distribution of material media making up an antenna or scattering structure 
by means of the constitutive equations 

D = e E + t H, B = vE + jx-H. (1 .5-1) 

Here the tensors € and /x have their usual roles; r and v allow for the description of possible mag- 
netoelectric properties of the medium [35]. 5 (This last bit of generality may provide future benefits 
and does not appreciably complicate the discussion.) The tensor parameters will of course in gen- 
eral depend upon position within the region of the antenna or scatterer considered; outside this 
region the set of parameters must reduce nominally to vacuum values: €= €o, /t = /x , T = v=0. 

In addition to a "given" or "original" system, described by the above equations, we must 
consider the adjoint system, which is described by the constitutive equations for the adjoint 
medium, 6 

D = €-E-*H, B = -tE + jaH, (1.5-2) 

where the superposed tilde denotes the transpose, and the tensors €, t, fJL, and v are those of 
the original system. As the equations show, "adjointness" is a mutual relationship: the adjoint of 
the adjoint system is the original system. 

Since the tensor parameters are essentially arbitrary, we should realize that the "original" 
system is in no way a theoretically preferred system; the designation is arbitrary but useful. 

The concept that a medium may be lossy, lossless, or even "gainy" is familiar. In an in- 
homogeneous medium these "dissipative properties" will in general change from point to point. 
In appendix B dissipative properties are more precisely defined and it is shown that these prop- 
erties are point-wise identical for mutually adjoint media. In view of the usual connection between 
dissipative properties and the concept of "realizability," we may say that mutually adjoint media 
are equally realizable. 

Nonreciprocal antennas are most commonly (if not invariably) so because of the use of fer- 
rites subjected to a static magnetic biasing field. In such cases the adjoint antenna can in prin- 
ciple be produced by reversal of the bias field. (We say "in principle" because in general no pro- 
vision is made for conveniently or precisely accomplishing the bias field reversal.) 

The scattering matrix for the adjoint antenna may be and is assumed to be defined with the 
same basis fields, the same reference surfaces, and altogether in the same way as that for the 
original antenna. Then, as shown in appendix A, the following generalized or adjoint reciprocity 
relations hold between the characteristics of mutually adjoint antennas: For the antenna-feed re- 
flection coefficients, 

Soo = S«; (1.5-3) 



r ' In some of the recent literature, media described by equations of the form of (1) are called "bianisotropic." However, optically active media, which have been 
studied for many years [36], satisfy equations of the above form with scalar parameters and are not anisotropic at all. Optically active media incidentally are reciprocal, 
and thus could form part of a reciprocal antenna or other device. 

6 Media related by (1) and (2) are called "complementary" media by Kong and Cheng [37]. The adjective "adjoint" seems more appropriate both nontechnically 
and technically: Maxwell's equations for the adjoint system can be written as the mathematical adjoint of Maxwell's equations for the original system. See appendix A. 
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for the transmitting and the receiving characteristics, 

TfoS^Cm, K) = -r ]m (K)S <]{) (m,-K), (g=l,2) (1.5-4a) 

T/o.So,(m,K) =-7i m {K)Si> (m,-K); (q= 1,2) (1.5-4b) 

and for the scattering characteristics, 

V m(K)S pq (m,K;n 9 L)=7}n (L)S<' )q (n,-L; m,-K). (p= 1,2; g= 1, 2) (1.5-5) 



The superscript "a" distinguishes quantities associated with the adjoint antenna. We observe 
that all the equations hold with S and S a interchanged. 

Scattering-matrix elements, such as Sio(wi,K) and Sn (ra,K,rc,L) , characterize processes 
associated with a pair of directions: 7 the direction of an incident wave and a direction of "ob- 
servation. " For the functions mentioned, the direction-pairs are no, e* and e/, e&, respectively. 
Reciprocity relates processes associated with two pairs of directions, obtained by reversing and 
interchanging the direction of incidence and the direction of observation. See figures 4 and 5. 

If the constitutive tensors obey the symmetry relations 



€, /x = fJL, and t — - 



(1.5-6) 



then, as may be seen directly from (1) and (2), the adjoint antenna and the original antenna are 
identical. In this case, if we use conventional terminology, we say that the original antenna is 
reciprocal; in the present context, a term such as self-reciprocal or self-adjoint would be less liable 
to ambiguity. In the self-adjoint case the superscript "a" is without force and may be eliminated; 
(3), (4), and (5) become expressions of properties of one and the same antenna (and reduce to 
results given with p = q = 1 in [1]). 

Equations (3), (4), and (5) are our basic expressions of reciprocity and adjoint reciprocity. 
As mentioned in the Introduction, these relations have found substantial application in research 
establishing the foundations of the extrapolation technique [6,7]. More immediate consequences 
are brought out in the next subsection and in section 4. 





Sb» 



s°,(-!<> ' 



FIGURE 4. Reciprocity direction-diagram for Sm and S({, 



7 This is in marked contrast to what is involved when spherical or other nonplanar waves are used to represent the fields in space. 
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FIGURE 5. Direction-diagrams for the scattering reciprocity relations. 
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1.6. Power Gain, Effective Area, and Polarization Indices 

Our primary purpose in this subsection is to exhibit some of the content and applicability of 
the PWSM formalism by defining some of the more familiar and conventional quantities of antenna 
theory in terms of the antenna scattering matrix quantities. In particular, we shall define the power- 
gain function G q and the effective-area function a Q associated respectively with the transmitting 
and receiving characteristics S 9 o and Sq q . We shall also define polarization indices p q o and p {)Q 
associated respectively with S^o and So q . All the quantities to be denned are of course functions 
of direction and all are meaningful and applicable in dealing with the coupling of widely separated 
antennas, the context of conventional antenna theory. It will be seen, however, that although these 
quantities are uniquely defined (in a given coordinate system) by the antenna scattering-matrix 
quantities, the converse is not true. Hence they are not adequate characterizations of transmitting 
and receiving properties for the type of theory of essential interest in this paper (and, of course, 
they say nothing about scattering properties of an antenna). 

Our definitions of power gain and effective area are consistent with the essential content of 
the corresponding IEEE Standard definitions. However, our definitions of receiving and of trans- 
mitting properties are formed wholly independently of each other — in contrast to a tactic used 
to some extent in the IEEE Standards. Reciprocity relations then appear only later in their proper 
role as theorems. 

The key relation appropriate for the present consideration of transmitting characteristics 
is the asymptotic relation 

E«(r) ~ - iysqo (RA/rJooe^/r. (1.6-1) 



Here 

s q o(K)a = b,(K), (1.6-2) 

where b q is the complete radiated spectrum in the sense of (1.2- 12a) and s^o is the corresponding 
"complete" transmitting characteristic given in terms of components of Sog by 

s,o(K) = [(-)^- 1 / 1 /y]S,o(KK)e„ +S^(2,K)ei . (1.6-3) 

The following relations are noted: 

k± -s„o(K)=0, 8g0 (0) =S,o(0), 

and Sgo(K) is the projection of s </() (K) on the k x , k : y plane. 

From (1) one easily finds for the power radiated per unit solid angle at large distances 

Pq=2 Fdy 2 |8«o(K)ao| 2 . (1.6-4) 

The Ko appearing here is the value of the admittance (eo/Mo) 1/2 for the ambient medium. The func- 
tion p q is a "power pattern" for the considered antenna; generally any function proportional to 
p Q , whether or not the factor of proportionality is known, is called a power pattern. 
The power gain (function) of an antenna is denned by 

r;,(K)=477/vn, (i.6-5) 
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where Po is the net input power to the antenna. It follows from (1.1-4) and (5) that 

g„(k) = *y ffi.y . d-6-6) 

i7o( i — poor) 

It should be noted that the power gain is a characteristic of the antenna under consideration, 
independent of the source used to excite the antenna. This means that the value of G (7 (K) is 
independent of the insertion or adjustment of a lossless tuner in the feed waveguide, whether or 
not this tuner is counted as part of the antenna (cf. discussion of eq (2.6-13) in [30]). 
Using the components shown in (3), we define the polarization index 

^ (K) -S Q0 (l,K)(-)^/y (L6 " 7) 

associated with the transmitting characteristics of an antenna. This spectral polarization index 
is definable more physically and more conventionally in terms of the components of the corre- 
sponding asymptotic E; by (1), it is just the ratio of the _L and 1 1 (or cj> and 0) components of this 
E in the direction of observation. Polarization characteristics are conveniently and fully described 
by the single complex number p^o; other polarization parameters, such as axial ratio and orienta- 
tion of the associated polarization ellipse, can be determined, if desired, from p^o. 

For an antenna functioning in a receiving mode, the counterpart of the power gain is the 
effective area or effective receiving cross-section, o~ f/ (K). Like the gain, this quantity is a scalar 
function of direction and involves a far-field concept — in this case that of an incident plane wave. 
It is here defined by 

fVmax=<r ff (K)Spoy, (1.6-8) 

where Pa, max is the available power at the antenna terminal and Sp oy is the magnitude of the Poyn- 
ting's vector of a plane wave arriving at the antenna from a given direction and providing a polari- 
zation match to the antenna receiving characteristic. The derivation of an expression for a q in- 
volves a number of intermediate results that are at least as important as <r q itself. 

Let the electric field of the wave incident on side q of the antenna be E-a exp (ik r)/27r; 
then from (1.3— 2a) or (1.3-4a), the emergent wave-amplitude at the antenna terminal is 

bo = Sooa + So g (K) • A. (1.6-9) 

(This expression differs from (1.3-8) only in that here we have not assumed ao= 0.) Thus, from 
(9), we see that the antenna, as excited by the incident spatial wave, appears at So as a source 
having a reflection coefficient Soo and a generated wave 6G = Sog(K) -A. By ordinary (microwave) 
circuit calculation the corresponding received power is found to be 

l (i-|r L |»)|So«(K)-^l 2 

^rec- o V0 pj L Q , 2 , (1.6-10) 

Z |1— J />3oo| 

where T L is the reflection coefficient of the passive termination at So. By setting r^ = Soo, we get 
for the available power 

Pa -~2 v " 1-ISool* (L6_11) 

To aid in the consideration of the polarization-related parts of the problem, we display the 
spectral vector 

A = oc\K\+ a 2 K 2 
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used here and originally defined in (1.2-18); we introduce the complete spectral vector for the inci- 
dent wave in the form 

a= [(-)%/y]aie||-fo:2e ± ; (1.6-12) 

and we define the polarization index 

for the incident wave. 

Further, it is convenient to define a receiving characteristic, so 9 , complementary to the com- 
plete transmitting characteristic s^o. The single essential requirement is that 

s 0(7 (K) • a= So 9 (K) • A (1.6-14a) 

be an identity in A. This leaves a possible ey-component of sor/ undefined; we are in fact free to 
require 

k^ -so 7 (K) =0. (1.6-1 4b) 

(Here, as before, the upper sign goes with 0=1, the lower, with q=2.) Equations (14) and (12) 
imply 

s 0a (K) = [(-)*y/*]Sog(l, K)e„ +So 9 (2, K)e ± . (1.6-15) 

One should note particularly that y/k appears here whereas k/y appears in the expression for s, /t) ; 
so,/ is in no way a "complete" vector of which Sog is a part. However, the relation 

8o 9 (0)=So 9 (0) 
does hold. 

Using the components shown in (15), we define a polarization index for the receiving charac- 
teristics 

Sog(2,K) 

'°' (K) - So„(l,K)(-)*y/* (L6 " 16) 

This parameter relates to the properties of a passive material structure; it does not directly char- 
acterize the elliptical path of a time-varying vector. It is a ratio of receptivities to components of 
polarization in an incident plane wave under the specified conditions. 

Observing that |a| 2 = 87T 2 Sp ylYo, we can now combine (11), (14a), (13), and (16) to obtain 



[1+jggPoJ 



(i+K| 2 )(i+ |p 0g | 2 ) 



i7T 2 7] {) \s oq \ 2 Spoy 

Fo(l-|Soo| 2 ) ( } 



The quantity in brackets is a "polarization mismatch" factor, which, according to the Schwarz 
inequality for complex vectors, attains its maximum value of unity for 



«i,(K) = po,(K). (1.6-18) 

Thus the condition for polarization match is expressed as a conjugate match of polarization indices. 
When (18) holds, we find from (17) and (8) 

47r^oiso^(K)_|^ n A 1Q v 

or q (K) — -T7-J- ■ ' (1.6-19) 

/o(l — |ooo| ) 

as the desired expression for the effective area. 
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The minimum value of Pa, incidentally, is zero: for the given antenna and for any given direc- 
tion of the incident wave, there is always a wave-polarization, w q = — llpo q , that is not received at 
all. The polarizations best received and not received are mutually orthogonal (in the power or 
Hermitian sense). 

The complete transmitting characteristic s^o and the complementary receiving characteristic 
So</ for a reciprocal antenna satisfy the reciprocity relation 

riokso q (K) = Yoys q o(-K), (1.6-20a) 

which follows from the basic relations (1.5-4). (In applying (1.5-4) it may be convenient to use the 
forms 771 = Yokly and 772 = Y y/k for the 77's, which were originally defined in (1.2-5, -6).) For mu- 
tually adjoint antennas, the relations corresponding to (20a) are 

rjoAsg^K) = Yoysqoi- K) , VWK) = F ys« (- K) . (1.6-20b) 

The power-gain and effective-area functions for a reciprocal antenna satisfy the well-known 
reciprocity relation 

(Tq(K) = ^G g (-K). (1.6-21-a) 

For mutually adjoint antennas the corresponding relations are 

o?(K)=^G,(-K), o-,(K) = £g«(-K). (1.6-21b) 

The results in this set are conveniently found as corollaries of (20). Relations of the type (21b) 
were noted by Harrington and Villeneuve for antennas containing gyrotropic media [38J. 

The reciprocity constraint for the transmitting and receiving polarization characteristics of 
a reciprocal antenna reads 

pog(-K) =-pgo(K). (1.6-22a) 

The corresponding relations for mutually adjoint antennas are 

P o y-K)=-p,o(K), p ,(-K)=-pV(K), (1.6-22b) 

which follow as further consequences of (1.5-4) or (20). 

An interesting corollary of the polarization matching and reciprocity theorems is that if the 
radiation from a reciprocal antenna in a certain direction is circularly polarized, the wave best 
received from that direction is circularly polarized in the same screw sense. A verification of this, 
though basically simple, involves a few key elements. The polarization-index definitions (7) and 
(13) are so fashioned that p r/ o(K) = i and w q (K) = i both represent right-handed circular polariza- 
tions 8 (for (7=1, 2 and for all K < /,). Here the behavior, and in particular the parity, of the unit 
vectors e M and e as functions of k is directly involved (refer to table 1 and especially to fig. 6). 
Note that e M , e 1? and e/ 4 - form a right-handed system congruent to e#, e y , and e z . The stated corollary 
follows upon setting p Q o(K-) =± i and applying (22a) and (18). 



8 See, e.g., Beckmann [39] or Hollis et at. [40]. For optical terminology (which differs) see. e.g.. Born and Wolf [33, p. 27|. 
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FIGURE 6. Unit vectors in and perpendicular to the plane of e z and kfor ± k, k real. 

2. Scattering-Matrix Analysis of Coupled Antennas; General Solution for System 

2 -Port 

We consider a system consisting of a pair of antenna systems operating in a homogeneous, 
isotropic, dissipationless medium, as shown in the highly schematic figure 7. We are primarily 
interested in this system as a transmission system, with one antenna transmitting and the other 
receiving. The complete treatment of a transmission system must include effects of scattering 
by both antennas, and thus automatically includes treatment of reflection systems, in which one 
antenna functions in both transmitting and receiving modes and the other antenna represents an 
arbitrary passive (linear) scattering object. 

For the description of the antenna on the left in figure 7, we apply (1.3-12). In the problem of 
interest there are no waves incident from the left (a 2 =0); the spectrum of waves going to the 
left, 62, is not involved in the process of solving the problem (but it is obtainable as a part of the 
solution). Thus (1.3-12) reduces essentially to 




FIGURE 7. Transmission system — schematic . 



bo = Sooao + Soiciu 



(2-1) 



These equations are set up with reference to terminal surfaces So and Si, the latter being at z=0 
in the coordinate system Oxyz. For the description of the antenna on the right, we again apply 
(1.3-12), using primes to distinguish quantities associated with this antenna. In the problem of 
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interest there are no waves incident on the right side of this antenna (a[ = 0) and the spectrum of 
waves going to the right, h[, is obtainable as a part of the solution. The needed scattering equa- 
tions thus are 

(2-2) 

K = -§20 a O+ S2 2 «2- 

For these equations the terminal surfaces are So and Si, the latter being z=d in the coordinate 
system Oxyz. 

The separation of the phase reference surfaces for the two antennas by the distance d implies 

a 2 {m, K) = 6i(m,K)c^ d , a x {m, K) = b' 2 (m, K)e iYd . (2-3) 

Now, with respect to the transmission path as an element of the system, the set of incident waves 
is represented by b\ and bl and the set of emergent waves by a\ and a! z . Thus from (3), the matrix 
description of this element is 

:;)=(m) 

where the elements of T are 7\K; L) = 8(k x — l x )d(k }l — I y )e iy(h)d . This is equivalent to the two 
separate transformations 

a { = fk, a 2 = fbu (2-5) 

We are now in a position to obtain a complete formal solution for the behavior of the trans- 
mission system under consideration. That is, we can obtain expressions for both bo and 6o, valid at 
arbitrary distances and including the effects of multiple reflections. (We can also formally determine 
the field in the transmission path.) We first consider transmission from left to right, assuming that 
the receiving antenna is terminated with a passive, reflectionless load. Using (2) and (5), we find 



oi=7 , &i=7S / 22a , 2 = 7S; 2 7 , 6i. (2-6) 

The operator TS 2 >T appearing here is the description of the receiving system, as a passive scatter- 
ing object, referred to the reference plane S\ of the transmitting antenna. Since this operator 
recurs frequently, we denote it briefly by R' . Substituting (6) in (1), we obtain 

6i = Sioa +Sn«'Si, (2-7) 

which (at least when written out more fully) is seen to be an integral equation determining b\. 
(It may be identified as an inhomogeneous, linear integral equation of the second kind.) The solution 
may be indicated formally by writing 

Si=(l-Sni?')- 1 Sioao. (2-8) 

This gives us the spectrum of outgoing waves in the transmission path; it includes both the simple 
plane waves and the evanescent waves. (a\ is now determined by (6); Ei/(r) and Hu(r) are deter- 
mined by (1.2-8).) We may obtain a more explicit but still formal solution to the basic integral 
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equation by the Liouville-Neumann method of successive substitutions. This leads to a representa- 
tion of the inverse operator in (8) in a series of iterated operators, 9 so that 

£i=[l+tf,i/£+(£i,l£) I + . . .]S, a . (2-9) 

The special virtue of this form is that the successive terms in the series correspond to successive 
round-trip multiple reflections between the transmitting and the receiving antennas. Of course, 
(9) is meaningful as an infinite series only if it converges in some useful sense. The domain of 
convergence will depend upon the "smallness" of the product SnR' , and it is worth noting that 
this product depends upon both Si i and S 22 (as well as upon the distance between the transducers). 

We complete this analysis by calculating the scattering matrix of the "system 2-port," which 
has its terminals at So and So and is defined by the equations 

b = Mnao+M l2 a' , 

(2-10) 
bo = M 21 a + M 22 a'o, 

(The properties of a transmission system are often conveniently embodied in this form.) Inasmuch 
as we have made ao = 0, solving for bo/ao and for b'o/ao yields directly 

Mi^Soo + Soi&'a-Siifl')- 1 ^, (2-11) 

M 2 i = S' m t(l-Siih')- l Sio. (2-12) 

A similar alternative solution with ao = and a' Q ^ yields 

m 22 =s; ){) +s; )2 r(\ -s 2 ' 2 R)- i s 2 ' (2-13) 

M l2 = S m t(\ -S; 2 R)-*S: {U (2-14) 

where R = TS\\T. These formulas were first given in [41]; formally identical expressions are ob- 
tained in the electroacoustics case [3|. Their general significance was mentioned in the Introduc- 
tion. Complete analytical solutions can be found in a highly specialized and idealized (but never- 
theless interesting) class of problems; in one such solution the associated Liouville-Neumann series 
is incidentally found to converge or diverge according as kd is greater or less than 0.87993310 . . . 
In the measurement technique to be described in subsection 3.1 (but not that in 3.2) we assume 
that the effects of reflections between antennas have been minimized and may be neglected. When 
such reflections are omitted, (12) and (13) become M 2 i = S' {)2 TSu) and M 22 = S' 00 , respectively. If 
the (passive) termination on the receiving antenna has reflection coefficient Ta, we obtain from (10) 

&o = F'S { y 2 TSioao, 

where F' = (1 — I /S 0() ) _1 . More explicitly, we have 

b' = F'ao J 2 S' 02 (m, K)S,«(m, K)e'"^/Iv. (2-15) 

(Equation (15) is essentially a basic and simple case of eq (43) or (46) in [1 J.) The integral appearing 
in this equation is called the transmission integral, and the scalar product in the integrand, 



9 The operator expansion in (9) is analogous to the finite-dimensional matrix expansion (1 —A) - l = 1 + A 4- A 2 + A 3 + . . . , which is valid if all the eigenvalues 
of A are less than unity in magnitude. 
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2 S' 02 (m,K)S 1 ,(m, K) ^Si 2 (K) • S I0 (K), ( 2 -16) 

m 

is called the coupling product. This marks the emergence of the central quantities involved in the 
antenna measurement techniques described here. 

In order to establish some of the content of (15) (but not for present applications) we note that 
the well-known Friis transmission formula can be derived from the asymptotic form of the equation. 
The asymptotic form in question is 

bo ~-2mkF'So2(0) - Sio(Q)a e ikd ld. (2-17) 

(This is analytically a version of (1.2~16a) or (1.2 — 16b) evaluated on-axis.) If we now calculate the 
ratio of the available power at the receiving antenna terminals to the net power input at the trans- 
mitting antenna terminals and use the Schwarz inequality, we obtain 

P„ 4*rd* ' (2_18) 

where G\ (0) pertains to the transmitting antenna, cr^O) pertains to the receiving antenna, and we 
have used (1.6-6, —19). Equality in (18) holds for polarization match. 

We mention one more important result, contained in (11). The first iterated integral in the Liou- 
ville-Neumann series for M n is the reflection integral. 



*(<*)' 



: f dKeMK)d$ 01 (K) • f S£ 2 (K, L) • S 1() (L)e^'>'r/L. (2-19) 



This is the simplest form of integral involving a scattering process. It can be interpreted as a mono- 
static radar equation, which, apart from multiple reflections between target and transceiver, is 
valid at arbitrary distances. Examples of the use of this equation may be found in [1]. In particular, 
upon setting S' 22 (K, L) = — 15(K — L), where 1 denotes the transverse unit dyadic, one obtains 
the interferometer reflection-integral of reference f 1 1. 

3. Determination of Coupling-Product Values 

3.1. Deconvolution of Transverse Scanning Data; Application of Sampling Theorem 

We can now quite easily give the analytical basis for determining coupling-product values from 
transmission data taken in a transverse plane. Let the required relative transverse displacement 
between the transmitting and the receiving antennas be denoted by a transverse displacement 
P = xejr+ ye y of the receiving antenna (fig. 8). By considering the phase kr (r=P + de z ) of the 
waves in the spectrum incident on the reference plane S' 2 of the receiving antenna, we see that the 
process of displacement simply introduces the factor exp (iK-P) in the integrand of (2-15). The 
expression for the received signal becomes 



aoF'j 



6i(P) = aoF' e ' K - p S' 02 (K) -S l0 (K)e*"dK. (3.1-1) 

Note that the quantity b' (P) is what is observed in the measurement process; it may or may not 
be simply related to E at the point (P, d). Inasmuch as (1) represents a Fourier integral trans- 
formation, its inversion is immediate: We write 

S,'«(K) -S,o(K)=D(K), (3.1-2) 

where D(K) is here an abbreviation for the determinate function of K given by 
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(3.1-3) 



The inversion of (1) is appropriately termed deconvolution (each of the factors S02 and S10 can be 
interpreted as Fourier transforms of certain physical fields). The term distinguishes the inversion 
of (1) from the inversion of the simpler equation (1.3-6) representing the Fourier-transform defini- 
tion of S 9 o. 




FIGURE H. Arrangement for transverse scanning: d is fixed, V is variable. 

Flic right-hand side of (3) is in fact determined (up to a phase factor) by the following measurable 
quantities: the distance d between the reference planes Si and S',; the reflection coefficients T/. 
and S'oo; and the magnitude and relative phase of b' {) (P)/ao as a function of P. (For convenience in 
measurement b' () (P) may be normalized in two stages, indicated by the expression 

/>;,(P) = fr;,(Po) fc?,(P) 

(U) do ^(P,))' 

where P is a selected fixed point.) 

An instructive hypothetical case in which our analysis correctly shows that no probe correc- 
tion, other than a known constant multiplier, would be needed occurs when the receiving antenna 
is considered to be an ideal electric-field probe. Such a probe can be described and treated ana- 
lytically as an elementary electric dipole antenna, assumed lossless and reciprocal. For such a 
probe it can be shown that the receiving response is b' {) (P) =CE(P) *e p , where e t , is a unit vector 
giving the orientation of the dipole and the absolute value of the constant C can be calculated. 
Evidently the assumed use of an ideal probe oriented in the x-direction, say, in (3), would effec- 
tively reduce (3) to the e ./--component of (1.3-6). 

As (1) shows, planar scanning (in the absence of multiple reflections) can be interpreted rigor- 
ously as a spatially invariant, linear, filtering process in the two-dimensional wavenumber domain 
(for filtering concepts see e.g., [53|). The concept of the action of a receiving antenna as a filter 
has been noted, e.g., by Brown [19], apparently with non-planar scanning in mind. The processes 
are substantially different in the planar and the non-planar cases, as becomes clearly evident 
upon consideration of the action of a highly directive receiving antenna. Nevertheless, the ability 
to account fully for the probe characteristics in the planar case suggests the use of a well-charac- 
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terized "large" probe to reduce data taking and processing effort. Measurements at NBS using 
such a probe have produced the expected beneficial results. The filtering process has also been 
studied by Joy and Paris [10J using digital filtering (after data taking) to simulate probe 
characteristics. 

Clearly the utility of the result (3) depends upon one's ability to evaluate the transform of the 
empirically observed b' (P) . Both least-square fitting and a two-dimensional form of the sampling 
theorem have been successfully used to evaluate the required transform of 6' (P) from data taken 
at the points of a rectangular lattice in the measurement plane [2, 8]. The application of the sampling 
theorem has become the method of choice, mainly because of the greater ease of computation, 
and will be described very briefly here. 

The essential requirement of the sampling theorem is that the function to be sampled be 
representable as the Fourier transform of a band-limited function. That />o(P) virtually fulfills this 
condition may be seen from (1): a band-limit Kb somewhat greater than A, and a distance (L may be 
chosen so that evanescent waves for all K> Kb are assuredly virtually zero in the measurement 
plane (e.g., with Kb = 1.05k and d=15\, attenuation at the band limit is approximately 260 dB). 
Bandlimiting within a smaller spectral region may result from the behavior of the product So2(K) • 
Sio(K) in individual cases (an early study is reported in [8]). If we assign band limits k\ =±2ttI\\ 
and k> — ± 27r/X 2 for k x and k y , respectively, a straightforward generalization of the usual one- 
dimensional theory, given in appendix C, leads to 

e -iyd 

D(K)= ^ «(Pn)e-*-«* (3.1-4) 

The vectors P/.s^irAiej+isA^e^with r,s= . . . — 1,0,1,2,. . .) define the measurement lattice, 
the quantities bl)(P rs ) are the (complex) values of probe output directly observed at the points of 
the lattice, and the summation goes over the points of the lattice. According to the sampling 
theorem, (4) is mathematically exact; that is, if the data [the bo(P r s)] were complete and exact, 
the result would be exact. (Although (4) is exact, it is not the "best possible" result; more advanced 
theory [42] shows that a rhomboidal lattice would be somewhat more efficient than the rectangular 
lattice.) The theorem requires an infinite sum, but in the applications thus far we have found that 
not even all values measurable above noise are needed. 

An important feature of (4) is that the highly efficient algorithm known as the "fast Fourier 
transform" is rigorously applicable to evaluate the sum. 

3.2. Received Signal as a Function of Distance and the Extrapolation Technique 

The extrapolation technique was introduced by Wacker and Bowman, has been described 
very briefly by Newell and Kerns [9], and more fully, in an experimentally oriented paper, by Newell 
et al. [11]. The theory and numerical techniques were developed by Wacker [6]. The papers men- 
tioned should be consulted for details; in the following paragraphs we give only a brief account 
of the main ideas and equations involved. 

The extrapolation technique requires that one observe b'o/ao as a function of antenna separation 
distance d, which is precisely defined by the choice of reference surfaces Si and S-> associated 
with the respective antennas (fig. 7). 

From (2-10) for the system 2-port, we obtain the expression 

b'o=aoM 21 l(l-M 22 T L ), (3.2-1) 

which is a precise and complete version of (2-15). Expressions for the elements M 2 \ and M 2 2 are 
given in (2-12, —13), from which, by a rather lengthy process, one finds for (1) as a function of 
d a series representation of the form 
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/'- //' V exp [i(2p+\)kd] ^ 

fto-ao/ 5 2," ,; 2/ , + i "2 A n<* "• (3.2-2) 

p = a q = 

We observe that the subseries of terms with a given p can be interpreted as the contribution of 
energy which has experienced 2p reflections or made 2/;+ 1 transits between antennas. In particu- 
lar, the subseries with p = involves no reflections and is the expansion of the transmission integral 

(2-15): 

¥(d) = £ V (doo + Aoid- 1 +A 2d-z+. . .). (3.2-3) 

a 

It is of considerable analytical interest that this series is not merely asymptotic but actually con- 
vergent for sufficiently large d, under the main hypothesis that the two antennas involved be of 
finite size. If d is measured between centers of spheres, of radii r and r\ each circumscribing one 
of the antennas, then d > do = r+r' is sufficient. A generally sharper, but more complicated, pre- 
scription for do can be given [6, 7|. 

By comparison of (3) and (2-17) we see that 

Aoo = -2mkS Q2 (0) • S,„(0). (3.2-4) 

Hence determination of the leading coefficient in (2) is tantamount to the determination of the (on- 
axis) value of the spectral coupling product. The basic idea of what we may call the conventional 
measurement method is simply to have d large enough to make other terms negligible compared 
to the leading term of the series. The basic idea of the extrapolation technique is to observe /> as 
a function of d and to fit this function with as many terms of (2) as may be significant, and so to de- 
termine a good value for //<)<> in particular. This enables one to cope with proximity effects and with 
multiple reflections between antennas. 

4. Utilization of Coupling-Product Data 

4.1. One Unknown Antenna (Transmitting or Receiving) 

In this subsection we provide basic equations for the use of coupling-product data in a basic 
antenna measurement situation: the measurement of an unknown antenna (transmitting or re- 
ceiving) with the requisite known antennas or antenna. Some of the concepts and notation estab- 
lished will be used in the next two subsections (which provide a partial answer to the question, 
"How does one obtain the first known antenna?"). 

In this subsection "nothing more" than the algebra of two (complex) equations in two (complex) 
unknowns is involved. The only likely case involving non-uniqueness and compatibility conditions 
is discussed. (In the following subsections the linear equations have to be reinterpreted as 
quadratics.) 

It will suffice to consider only the case in which the unknown antenna is transmitting; the dis- 
cussion of the other case would, of course, be analytically very similar. 

We do not make simplifying a priori assumptions concerning symmetry or polarization (or 
other) characteristics of the antenna to be measured. Quite generally, then, we require measure- 
ments to be made with (at least in effect) two receiving antennas, A and B, having suitable known 
receiving characteristics, A .>(K) and Bo2(K). From such measurements, the values of the coup- 
ling products 

A 02 (K) -S 1( )(Iv)=Z),(K), 

(4.1-1) 
B 02 (K) -Sio(K)=ZMK) 
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are to be determined for the desired values of K. For each chosen, fixed value of the parameter K, 
we have two (complex) equations for the two (complex) components of Sio(K). 

Complete solvability of (1) requires that the vectors A02 and B02 be linearly independent at the 
value of K considered. A measure of the linear independence of the two vectors is given by the 
following expression for the normalized squared-magnitude of the determinant of (1): 

|A|2 ' |A -« 12 (4,-2) 



|A|«|B|« |A|»|B|» 

As the Schwarz inequality shows, this quantity ranges from the value zero when B02 is proportional 
to A02 to a maximum of unity when B02 * A02 vanishes. In other words, power orthogonality repre- 
sents the extreme case of linear independence. Examples are linear polarizations at right angles 
and left- and right-circular polarizations. 

An explicit, coordinate-free solution of (1) may be obtained with the use of the set of vectors 
reciprocal to A02 and B02. The reciprocal set a, /} is defined by 

a A ,,2= 1, a-Bo2=i3-Ao2 = 0, /2B«> 2 = 1 

(assuming the required linear independence), and 

Su> = D A a+D B p (4.1-3) 

as is easily verified. The algebra is summed up in the statement that D A and D B are the covariant 
components of S10 with respect to A02 and B02 as base vectors. 

In much of what follows the use of %, y components relative to the fixed basis e#, e^— rather 
than the 1, 2 components relative to the variable unit vectors k x , k 2 — is indicated. The considera- 
tion of geometric rotations and symmetries is appreciably complicated by the dependence of the 
k's upon K. Also, the use of e x , e y automatically takes care of the matter of defining the k's on- 
axis. The required coordinate transformations will be governed by the relationship of the unit 
vectors 

Kx = ce x + se y , K 2 = z — se x + cey, (4.1-4) 

where c = cos (/) =kJK and s = sin (p = kyjK. 
We introduce the abbreviations 

4r = >WK), B x =B* 2x {K), S, = Si te (K), (x = x,y) (4.1-5) 

and write (1) in component form 

Ajc^JC I AyOy = L)\ , 

(4.1-6) 

B X S X "+" BySy = D[{. 

Instead of using two intrinsically different antennas A and B, in many cases it may be possible 
and convenient to use one antenna in two orientations, differing by rotation around the z axis by 
90 degrees, say in the direction of x to y. If antenna A is so used and so rotated, we obtain for the 
equivalent of antenna B 

Bo2 X (kx<> ky) — — di)2y(ky, — h x ) 

(4.1-7) 
Bo2y{kx , ky) — A() 2x (ky , — kx) • 
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These equations express the rotation of the vector field A02 corresponding to the rotation of the 
antenna that it describes. 10 They lead to a modified version of (6), which we shall discuss in the 
particularly interesting case of evaluation on-axis (K = 0). If we let A x =//<>-., ■(()) , S x Siox(O), 
etc., the modified form of (6) is 

(4.1 8) 
-^0+^250 = ^(0). 

The determinant of the system (8) vanishes if and only if A° x =±iA Q i that is, the determinant 
vanishes if and only if the response characteristic A02 is "circularly polarized" at the point K=0. 
This gives us a hint as to the special advantages of the use of circular polarization components for 
equations of the above form. We introduce circular polarization components for on-axis quantities 
in the following manner: For the transmitting characteristic 

S 10 (0)=So e + +S°_e- (4.1-9a) 

and for the receiving characteristic 

Ao2(0)=^ , je + +^ ( le_, (4.1-9b) 

where the superposed bar denotes the complex conjugate (as usual) and 

e f = (e x + ie y )lV2, e_ = (e^ — je y )/V2. 

(The ordered pair of vectors e f , e in (9b) is reciprocal to the ordered pair e f , e_ in (9a).) The 
formulas for transformation to circular components are thus determined as 

S°= (S° + +S°_)/V2, S° y =i(S%-S _)lV2, (4.l-l()a) 

A^=(A»+A»)/V2, A%=-i{A\-A'l)lV2. (4.1-10b) 

The scalar products in (8) become 

A<\S<> + + A°_S<L = Da(0), (4.1-1 la) 

iA° + S + -iA<LS<L = D B (Q). (4. 1-1 lb) 

As is apparent, we have chosen the notation so that A f and A_ represent receptivities to the corre- 
spondingly labeled circular components of S10. The coupling is of course consistent with the 
polarization matching theorem (1.6-18). Equations (11) show clearly that if the receiving antenna 
were to respond to only one circular component of polarization, then that component, but not the 
other, could still be measured using one or the other of the two equations. If both equations were 
to be used, the compatibility condition (required by the vanishing of the determinant) should be 
satisfied within experimental error. 

4.2. Generalized Two-Identical-Antenna Techniques 

The technique to be described is formulated for nonreciprocal antennas, assuming that one 
possesses the adjoint of the antenna to be measured. If the antenna to be measured is reciprocal 



10 Concepts, analytical tools, and notation for the rotation of vector and tensor fields are discussed on p. 272 of [43]. 
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(= self- adjoint), then the assumption is that one possesses duplicate antennas. 11 This is no doubt 
the more likely case but the more general formulation can be given with essentially no extra alge- 
braic complication. Under the main assumption we are assured that the receiving characteristic 
of one antenna will be related to the transmitting characteristic of the other by reciprocity relations 
whether or not the antennas are reciprocal. 

Granted the assumption of duplicate or mutually adjoint antennas, no additional assump- 
tions are required to permit formulation and solution of equations for on-axis gain and polarization 
characteristics of both antennas. Additional a priori information, ordinarily qualitative, is required 
only for resolution of square-root sign ambiguities. 

Certain commonly occurring types of symmetry permit one additionally to obtain solutions 
for off-axis values of gain and polarization. 

For definiteness we assume transmission from left to right (as usual); limit the discussion to 
the determination of the right-side characteristics of the two (in general distinct) antennas labeled 
S and S a , say; and choose to formulate equations for the direct determination of transmitting charac- 
teristics, leaving receiving characteristics to be determined by reciprocity. Under these ground 
rules, the remaining problem consists of at most two parts: 

(a) Transmit from S to S a ; formulate equations for Si , find Sgfj by reciprocity; 

(b) Transmit from S a to S; formulate equations for Sf , find Soi by reciprocity. 

If reciprocity (in the ordinary sense) applies, the superscript "a" is without effect and may be 
omitted; the two cases reduce to one. 

Inasmuch as the algebraic problem is in all cases substantially identical, it will be sufficient 
to consider only case (a) explicitly. 

Both antennas are initially to be described in the same orientation and position relative to 
the fixed coordinate system Qxyz. When one of the antennas is placed and oriented to serve in 
reception, its description relative to fixed coordinates will be changed accordingly. Indeed, the 
phase factor exp(iyd) introduced by the axial translation can be regarded either as a modification 
of the receiving characteristic of the receiving antenna or as a property of the transmission path. 
We make the latter point of view explicit by referring the description of the receiving antenna to 
the shifted coordinate system O'xyz', where 0' is at point (0,0, d) in the original system Oxyz 
(fig. 7). 

We shall need to consider the receiving antenna in two receiving orientations, differing by a 
90 degree rotation around the z-axis. The operative characteristics of the receiving antenna in 
these two orientations will be distinguished by single and double primes. 

Let the adjoint antenna be rotated into the first receiving orientation. This requires 180 degrees 
rotation around a transverse axis, say the y-axis. The operative receiving characteristic of the 
rotated antenna is then 

Sfc'(K)=P ir 2Sfi(K), (4.2-1) 

where P y2 is a notation for the transformation (of SgJ produced by the prescribed rotation. Now 
the transformed function Pj^Sfl is related by reciprocity, (1.5-4), to the similarly transformed 
function P^Sio. Thus, if the reciprocity relation is written in vector form, we have 

i?oSft'(K) = i|-P^Sio(-K); (4.2-2) 

and if it is written in x, y component form we have 



11 The idea of using a reflecting surface or mirror to produce an image antenna is not fully applicable, even if the antenna is reciprocal: Coupling-product data 
can be obtained by extrapolation or by conventional techniques, but not by transverse scanning. The coupling-product equations can be formulated and solved for 
on-axis values | provided on-axis polarization is known and not circular. It should be noted that the image antenna, being a mirror image, cannot in general be considered 
an identical antenna even with a perfectly reflecting surface of infinite area. 
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-Si te (*„-*,) \ 

(4.2-3) 

Swy (/,.,. -/.■„) / 

Here the square matrix is determined by (4.1-4) and the column matrix on the right contains the 
%, y components of Py2Sio("~K). We may now evaluate the coupling product Sio(K) • S ( ',(/(K) 
= D'(K). Some degree of abbreviation is indispensable; we use 

S x =Sio x (K), S y =S, y(K), S|=Si to (*x,-Ay), Sj=Si 0y (Ax,-Ay), (4.2-4) 

and obtain 

- (i7 1 c 2 + T7 2 5 2 )S iZ .5]:-f (171 — ij2)5cSxSg.H- 

-(T7i-T7 2 )5cS 1/ Si?+(7) 1 5 2 -hr)2C 2 )S 1/ S^=77o/>'(K). (4.2-5) 

This is the first of the desired "measurement equations" relating the mathematical expression 
of the coupling product to its empirically determined values. It is interesting that the expression 
is invariant with respect to the interchange of k y and — k y ; the empirical D'(K) should also have 
this symmetry. Further, since D' (K) is related to b' (P) by (3.1-1 ).the b' {F) data should have 
the corresponding property of invariance with respect to the interchange of y and — y. These 
general constraints should be experimentally useful. 

To obtain the second measurement equation, we rotate the receiving antenna, as described 
by (1), 90 degrees around the z-axis in the direction .v to v (cf. (4.1-7)). Using the notation P z a 
for this rotation and applying the reciprocity relation, as in (2), we obtain 

T ?0 S b a 2 "(K)=7 ? -P,4P y2 S,o(-K). (4.2-6) 

In x, y component form this is 






-Siox(k y ,k x )J 



(4.2-7) 



Using the abbreviations 

Sj? = Si te (A y , k x ), S d y =S m (k y , /,-,), (4.2-8) 

as well as S x and S y in (4), we find for the coupling product Sio(K) • $\y>{ k) = /)" (K) the expres- 
sion 

-(Vir 2 + V^' 2 )S,S<i-( Vl - V2 )scS J .Si + 

-(7 ?1 -T ?2 )5cS y S^-(T ?1 5 2 -fT) 2 C 2 )S y S|= 7]oD"(K). (4.2-9) 

In this case the coupling product is invariant with respect to the interchange of k x and k y , and 
again this constraint should be experimentally useful. 

In the remainder of this subsection we discuss briefly some conditions and methods for 
determining components of S U) from (5) and (9), assuming that the D's are given for the values 
of K of interest. When evaluated for K ¥= 0, the six quantities S x ,S y , S v x , S v y , S x , and S* appearing 
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in these two equations are in general distinct and unknown; obviously, some specialization or addi- 
tional data are required. We consider two cases: evaluation on-axis, and a simple type of symmetry. 
(i) Evaluation on axis — Evaluated on-axis, the six unknowns reduce to two 

S% = S 10 * (0) , Sg = Sioy(0) , (4.2-10) 

and the coupling-product equations reduce to 

- (Sn 2 + (Sg) 2 = D'(0)tioIYo, (4.2-lla) 

-2S%S Q y = D"(0) V olY o . (4.2-llb) 

(The key to this reduction is the observation that Yq is the common value of r?i and r]i on-axis.) 
Subcases under this case occur if the polarization on-axis is considered known. One or the other of 
(11) will suffice, no matter what that polarization may be. 

Suppose, for example, that p = S®/S® is considered known and not equal to ± 1; then from 
(11a) we may obtain 

^°^ ' (4.2-12) 



F„(p*-1) 

This, together with .S!) = pS%* gives us the on-axis pattern vector in terms of D' (0) (up to a sign). 
From the expression (1.6-6) for power gain we find 

(1 - |Soo| 2 )|p 2 - 1 1 

We can obtain an interesting form for this result by expressing D' (0) in terms of the integral 
of b[) (P), as in the deconvolution relation, (3.1-3). Thus, 

G, (0) = £ '/"'ffi 1 £il±-j ±- | / b' n (P) dP |. ,4.2-14) 

7t i — |£oo r IP — M l a <>| 

(This result, as well as the version of it with p = 0, has been presented previously [41, 44].) The es- 
sential simplicity of the result is somewhat obscured by the presence of the mismatch factors. 
If we assume a polarization match (p=0, ± /, or °°), a conjugate impedance match (r,=S ( ' )0 ), 
and Soo = S' {)i) (as appropriate for mutually adjoint or for identical antennas), we have 

Ci(0)=~~| f b' (P) dP\. 

77 |a | 

For the effective area, using (1.6 — 21b), we obtain the remarkably simple expression 

crf(0)= i -i r | fb' Q (P)dV\. (4.2-15) 

|a | 

Equation (13) has been successfully applied experimentally [8|, and (15) has been tested analytically 
in a special case. In the analytical test the two identical antennas were taken to be ^-oriented ele- 
mentary electric-dipole antennas, assumed lossless and reciprocal. In this case one does indeed 
obtain the expected result cri(0) = 3X 2 /(87r). 

More generally, the polarization is not known and it is necessary to solve (11) as simultaneous 
quadratics. A solution in terms of circular polarization components is convenient and useful. 
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Using the definitions in (4.1 -9a), one obtains 



(S° + ) 2 -■$■ [D'(0)-iD"(0) 



(S°J 2 = -■£- [D'{0)+iD"(0)]. 



(4.2-16) 



2F„ 
From these equations we obtain four pairs of values for the x, y components: 

S%=±± J^ ( VD'(O) + iD"(Q)±VD'(0)-iD"(0)). 

i— (4.2-17) 

S»,= ± 7j y-p 9 ( VD'(O) +iD"(0) + V£>'(0) - iD"(0) ). 

(The double signs are correlated vertically but not horizontally.) Which solution-pair pertains to 
a given measurement cannot be determined from the equations alone. It would seem that ordi- 
narily the overall plus-minus signs should be of no significance. The remaining sign-choice does 
affect the determination of the polarization index p = SyS {) r . If we write p (u) and p U) for the values 
associated with the upper and the lower signs, respectively, we find 

pM p U) = - |. (4.2-18) 

Witli the aid of this equation a modicum of a priori information about the magnitude and/or the 
phase of p should ordinarily be sufficient to resolve the ambiguity. For example, if the antenna is 
known to be approximately linearly polarized (on-axis) in a certain direction (we may choose the 
y axis in that direction), then \p\ is distinctly greater than unity and one would choose p (,/) or 
p (/) accordingly. However, if the polarization is nearly circular, the difference between the two 
indices becomes relatively small and the choice could be difficult. 

The squared magnitude \S%\ 2 + \S°y\ 2 , which determines the on-axis power gain, is unam- 
biguous. In fact, for this gain we find 

C|(0)= 2^M|^(0)-H^(0)| + |/>'(0)-»P"(0)|] 

1 — |ooo| 

(ii) Additional solutions permitted by symmetry. Certain types of symmetry enable one to de- 
termine certain off-axis values of the unknown functions. We consider one type of symmetry of 
frequent occurrence, exemplified by the fields of (a) rectangular waveguide open-ended or with a 
pyramidal horn, fed by the TEw mode in the waveguide; (b) circular waveguide open-ended or with 
a circularly symmetric horn, fed by the TE\\ mode in the waveguide; and (c) a transverse electric 
dipole. In these examples the symmetry may be analytically specified in terms of spectral com- 
ponents by 

S\0x\ A\ r , tCy) = SlOx(kxi ky) = O K).r ( K x ,k y ) , 

(4.2-20) 

S\0,f{ — k X ,ky) = Sl0y(k X , — ky) = S i<) y(k X ,k y) . 

(Here a specific orientation of structures and fields has of course been assumed.) We note in 
particular that S\o x must vanish on the coordinate axes. Hence, the four unknown functions in- 
volved in (5) are reduced to the single one, Sio y (k x ,0) on the line k y = 0, and to the single func- 
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tion Sio y (0,ky) on the line k. r = 0. From (5) we immediately obtain the separate equations de- 
termining these two functions, 

(4.2-21) 

7l 1 [SlOy(0 9 ky)y=7loD'(0 9 ky). 

Next we notice that on the diagonal line k y = k x , we have S^ = S x and S d = S y . Additionally, the 
symmetry furnishes S x = — S x and Sy = S l y. Hence, of the six unknown functions involved in (5) 
and (9) essentially only two distinct ones survive, 

u x = Siox(k x , k x ), u y = Sio y (kx, k x ) , 

say. Further, the two quadratics can easily be solved for these unknowns, as the following inter- 
mediate results show: 

1)2 (a, - u y y = rio[D'(kx, k y ) + D"(k x ,k x )] 9 

(4.2-22) 
17, (u x + u y y = no[D' (*„ k x ) - D"(k x , k x ) ] . 

Hence in this case we can obtain solutions on the four lines, t-0,^=0, and ky = ± k x . 

4.3. Generalized Three-Antenna Techniques 

In this subsection we discuss the analysis involved in techniques for determination of both 
power gain and polarization using three unknown (dissimilar) antennas. We require an antenna 7\ 
to be used only in transmitting; an antenna /?, to be used only in reception; and an antenna S to be 
used in both receiving and transmitting modes. We do not need to inquire whether either of the 
antennas T and R is reciprocal or even capable of operating in a "reversed" mode. We do require 
either that antenna S be reciprocal or, if not reciprocal, capable of being "switched" to become its 
own adjoint S a . 

It is interesting and important that some kind of reciprocity 12 is indispensably required a priori 
information. This requirement cannot be avoided by increasing the number of antennas involved — 
even to the extent of using all possible combinations of n transmitting and m receiving antennas. 
But when the reciprocity requirement is met, as with a reciprocal or switchable antenna S, 3 
antennas are sufficient. (Note that in the switchable case, S and its adjoint do not coexist. Actual 
simultaneous possession of both S and S a would distinctly change the character of the measurement 
problem: In this case one could use the generalized 2-identical-antenna technique of the preceding 
subsection.) 

Of course the use of 3 antennas in a measurement scheme is not in itself new. The use of 3 
antennas in the roles of T, R, and S (S being reciprocal) is recognizable in a discussion in Vol. 12 
of the MIT Radiation Laboratory Series — material originating in the early 1940's [45]. However, 
in that discussion simplifying assumptions regarding polarization were made and the analytical 
problem was reduced to the use of scalar equations involving gain only. Analogous schemes using 
three electroacoustic transducers coupled by a fluid medium are well known [46]. The application 
of the requisite electroacoustic reciprocity relations was introduced by Maclean [47] and by Cook 
[48] in 1940 and 1941, respectively. 

In this discussion we consider only on-axis values of antenna characteristics and refer to these 
quantities as the antenna characteristics. Determination of off-axis characteristics by three- 



12 Conceivably some linear relation, other than that provided by ordinary or adjoint reciprocity, between receiving and transmitting characteristics might he kn 
i priori. This would not in general lead to the relatively simple equations that are obtained under the assumption made. 
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antenna techniques has not yet been seriously considered, to our knowledge. (The corresponding 
electroacoustics problem, which is much simpler, has recently been completely solved under 
minimum assumptions in work as yet unpublished.) 

When antenna polarization characteristics are fully taken into account, as in the present 
discussion, three-antenna techniques for the determination of on-axis characteristics require 
less a priori information than any other technique. Thus versions of the three-antenna technique 
have been the method of choice in several critical applications [e.g., 9,11]. 

The problem now posed is taken to be the determination of the receiving characteristic 
Roi (0) of the antenna R, the transmitting characteristic Tio(0) of the antenna T, and the four 
* 4 right-side" characteristics Sio(0), Soi (0), S" (0), S^ (0) of the antennas S and S a . We further 
choose to eliminate the receiving characteristics of S and S" by means of the reciprocity relations, 
which for K = are simply 



T7oSoi(0) = KoS;' (0), ijoSg, (0) = K S,o(0). 



(4.3-1) 



In this discussion we use the notation 



R* = Ro2x(0), S, = .s,o.,.(0), s:; = S;v,.(0), r,= r,o*(0) (x = x,y) 



(4.3-2) 



for the components of the four remaining vector unknowns. Certainly, to determine these H 
(complex) quantities we need a system of 8 (complex) equations; the "core" of the problem, however, 
turns out to be the solution of six simultaneous quadratic equations for six of the unknowns. (This 
is the complete solution if antenna S is reciprocal.) The remainder of the problem requires only 
the solution of two linear equations in two unknowns; this part will be called the "supplementary" 
problem. 

The complete problem requires coupling-product data derived from transmission between 
antennas paired as shown in table 2. Transmission is from left to right (as already implied in the 
notation), and for each antenna-pair considered, the receiving antenna is used in two orientations 
differing by 90 degree rotation around the z axis. The first two columns in the table indicate two 
ways of utilizing the same data and are related by the interchange of the roles of S and S". If the 
antenna S is reciprocal (S = S a ), the first two columns coincide and reduce to the third; only the 
core problem remains. It will be sufficient to formulate and solve the algebraic problem represented 
in the first column of the table; the corresponding formulations and solutions may be obtained for 
the second column by interchanging S and S a and for the case of reciprocal S by eliminating the 
superscript "a" (and with it the equations thus rendered superfluous). 



Table 2. Antenna pairings 



For "core" problem: 

For "supplementary" problem: 



S = S" 




For the transmission from T to R, we have 

T,R, + TyRy = D' RT , - T X R y + T ' y R x = D" RT . 



(4.3-3) 



These equations are an instance of (4.1-8) (except that here characteristics of both antennas are 
unknown). For the transmission from S to R, we have 



SjR X + S yR y — D' RS , — Sj-Ry+ SyRjc — D" RS . 
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(4.3-4) 



These equations are similarly an instance of (4.1-8). For the transmission from T to S a , we have 

- TVS* + T y Sy= r, D' ST IY , - T,S U - T y S x = tjo^/Fq. (4.3-5) 

The receiving antenna, S a , is here treated in the same way as was the receiving antenna in the two- 
antenna technique described in the preceding subsection (cf. (4.2-10)). This completes the formula- 
tion of the "core" problem. The solution is conveniently accomplished with the aid of circular 
polarization components, which are defined as in (4.1-9). Equations (3), (4), and (5) transform by 
pair to 



T-R- — -(D' RT + iD" RT ) = S/jr, 



(4.3-6) 



(4.3-7) 



T+R + = -{D' RT -iD" RT ) -A* r ; 
S-R- = ±(D' RS + iDl s ) S J W , 

S+R + =\{D' RS -iD" RS ) 3 A M ; 

T + S + = -± V o(D' ST - iD" ST )IY {) - - tjoAWFo, 

(4.3-8) 
TS- = -±Vo(D' ST + W" ST )IYo - - thSstIYo. 

For the squares of the individual circular polarization components one finds without difficulty 
(but note that neither this set of equations nor the preceding sets can be generated by straight- 
forward advancement of subscripts) 

R2j= _Yo^hs at _ s _^ 0) 

(4.3-9) 
R%=-— -"'-" ; (T + S + #0) J 



\ (4.3-10) 

*0 ^RT 

n^-W-^T^As+R+rto) , 
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(4.3-11) 



Yo Zrs 



The inequalities shown in parentheses (which in effect prohibit circular polarizations) must be in 
force if indeterminacy in the respective associated equations is to be avoided. If circular polari- 
zations do occur, then obviously many special cases are possible. For example, if the transmitting 
characteristic of antenna S is (very nearly) right circularly polarized (S y = iS x ), then (very nearly) 
S- = 0, S/es=0, X.s7 = 0, and R- and T- become experimentally indeterminate. However, the 
determination of the other circular components of R and 7\ and the determination of the trans- 
mitting characteristics of antenna S itself, would present no special difficulty. In practice, if a par- 
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tial solution such as this is not sufficient, one might include additional antennas in the scheme. 
Antennas capable of being switched between right and left circular polarizations have been used 
for similar reasons, usually in less complicated circumstances. 

It remains to dispose of what we called the supplementary problem earlier in this discussion. 
If antenna S is reciprocal, there is no supplementary problem; if S is not reciprocal, it remains to 
determine S? (0). As suggested in the first column of table 2, this can be done by measurement 
with the receiving antenna R, for which calibration is presumptively available as a result of the core 
problem. Equations (4.1-1 1) are applicable; in notation adapted to the present context those equa- 
tions become 

#+S« + R-S a _= D' RS a 

iR + S a + -iR~S<L=Dl s «. 

Here complete solvability requires R + R- ¥" 0. 

Effective areas and power gains for antennas R, S, S a , and T(as pertinent) are readily expressed 
in terms of circular polarization components, which by reference to (4.1-9), are seen to be normal- 
ized so that 

l/M 2 + l^l 2 =M + l 2 +M-l 2 - 

Thus, for example, we obtain for antenna S from (10) and (1.6—6) 



ooo 



A/jt 



Z*RT 



(4.3-12) 



The corresponding effective area is given by (1.6— 21a) or (1.6— 21b), as appropriate. 

Algebraic values for the individual linear polarization components may be obtained from 
(9), (10), or (11), as desired. For example, for antenna S 



ZiRS^'ST . /A/f.s-A.s'y 



v 



ZiHT > Art 

(4.3-13) 



2J0 / /S/f.yS.sT — /A/f.yA.sT 

Yo\ V Zrt + V A rt , 



where the double signs are correlated vertically but not horizontally. As far as sign ambiguities 
are concerned, these equations have the same algebraic structure as (4.2-17) and the earlier dis- 
cussion is again applicable. 

5. Appendices 

5.1. Appendix A: Reciprocity Theorems 

The domain of the electromagnetic fields under consideration is the source-free region V 
bounded externally by the surfaces F\ and F% and internally by the closed surface So + S a , which 
encloses the source or detector associated with the antenna considered. (See sec. 1.1 and fig. 1.) A 
passive antenna is a scatterer (sometimes called a "loaded scatterer"); if the structure is merely 
a scatterer, the surface So + S a is irrelevant and may be disregarded. 

We shall write Maxwell's equations in a form especially suited to the purposes of the present 
discussion. We associate with the given system, described by the constitutive equations (1.5—1), 
the "Maxwellian" operator 
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\— ioji>- + Vx —ioj/ji- ) 
and place the vectors E and H in the column matrix 



-(h> 



so that Maxwell's equations for a possible field in the system are expressed by 

McD = 0. (A-l) 

Using the constitutive equations (1.5-2) we similarly obtain for the Maxwellian operator associated 
with the adjoint system 



M 



/ ia)€ • — icov • + V X \ 

V'oot*+VX —iojfJL' ) 



where M is in a natural sense the transpose of the operator M. Maxwell's equations for the fields 
in the adjoint system take the form 

McD = 0. (A-2) 

The adjoint system represents a possible physical embodiment of the mathematical concept 
of adjoint differential expressions, which may be defined as follows. Let 

*'=(£)' *" = (£)' < A ~ 3 ) 

where F\ G', F", and G" are arbitrary differentiate vector functions denned in V. Then we can 
associate a unique adjoint differential expression M" with M by requiring that 

<£>" -MO' - & • M«4>" = D, (A-4) 

where D is a divergence expression, be an identity in 4>' and 4>". This leads to 

M" = M (A-5a) 

and 

D=V-(F , X G"-F" xC). (A-5b) 

Thus the adjoint operator turns out to be the transpose. 13 

We may now easily obtain the basic theorem for our purposes. If we replace F\ G' and 
F", G" by E', H' and E", H", which satisfy (1) and (2), respectively, we have 

V-(E'x H"-E" xH') = (A-6) 

throughout the region V. This is a generalization of the well-known Lorentz relation. It may very 
easily be modified to include the effects of both electric and magnetic current sources embedded 
in V [37]. Here, however, we are interested in excitation of the systems by means of incident waves 



"To avoid possible confusion, we emphasize that the adjoint arising here is definable essentially by transposition (without complex conjugation); it is not the 
Hermitian adjoint that is frequently useful in other physical problems. For discussion of the theory of adjoint differential expressions see especially Lanczos [49] or 
(for ordinary differential equations) Courant-Hilbert [50]. Equation (5a) incidentally contains the nontrivial result that the operator V X is self-adjoint. 
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(in space or waveguide). With or without current-source terms, (6) characterizes what may properly 
be called adjoint reciprocity. 

If the constitutive tensors obey the symmetry relations 

€ = €. M = M- andT = -^ (A-7) 

then (as may be seen from (1.5-1 , -2)) the adjoint system and the original system are identical (and 
M a = M=M). In this case we say that the linear differential operators and the systems are self- 
adjoint— and ordinary reciprocity obtains. Equation 6 still holds; E\ H' and E", H" may be in- 
terpreted as distinct electromagnetic fields in one and the same system. 

To apply the generalized Lorentz relation, we first take the volume integral of the expression 
over the region V and use the divergence theorem. This yields 



f L. no dS=2 | 4-n„rfR. (A-8) 



where L=E'XH"-E"XH' and the unit normals are inward on So and outward on F q (as 
prescribed in subsections 1.1 and 1.2). Next we substitute the modal representations of the fields 
E\ H' and E", H" on So, F\ and F>* using single and double primes to distinguish the spectral 
variables associated with the respective fields. One obtains after some analysis the generalized 
reciprocity lemma 



r)o(a'obo — aj&o) 



£ [ £ [a£(rt,-L)6;,(/i, L)-ai(n, L)6J(n,-L)] t?„(L) rfL. (A-9) 

p •'I. n 
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From this lemma, by making suitable special choices of excitation, we shall obtain the desired set 
of relations between the elements of the scattering matrix of the original system and the corres- 
ponding elements of the scattering matrix of the adjoint system. 

For immediate reference we write down the scattering equations for the original system 



b' = Sooflo "+" X I 2 S°p( n ' ^) a' p (n, L) rfL, 

p •'L ti 

b' q (m, K) = S q o(m, K) a + X I 2 ^QP^ m -> ^; ^' ^) 0p(ra, L) dL. 

and those for the adjoint system 

K = S&flS + 2 f 2 Sg„(n, L) a>, L) ,/L, 



(A-10) 



l> J L n 



;'(m, K) =S|o(m, K) «,',' + ^ f S S &( m ' K ' "■ ' > "/>< l ' ] dL 

l> •'L n 



(A-ll) 



The superscript "a" distinguishes quantities characteristic of tin ystem. The above sets 

of equations are instances of (1.3-2); each set is defined relative to the same basis fields and refer- 
ence surfaces. 

To find the relationships between transmitting and receiving characteristics, we let E', H' 
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and E ", H " be the fields corresponding to excitation of the respective systems by the incident waves 
represented by the following set of spectral variables: 

a' = 1 , a' p (n, L) = 0, 

ao = 0, (i pin, L) =8 pq dnm 8(L — K). 



From (9) we find 

and from (10) and (11) 



-b'oT)o= bq(m, — K) r) m (K) 



b q (m,K)=S q0 (m, K), b; = SV«, K). 

Hence [observing that r) m (—K) = r) m (K)] we obtain 

i)oS&(/», K)=-T?,,,(/v)S y() (w,-K). (A-12a) 

Similarly, by interchanging the patterns of excitation (or by interchanging the designations of the 
original and the adjoint system), we obtain 

r)oSo q (m, K)=-r) m (K)S q \(m,-K). (A-12b) 

To obtain the scattering reciprocity relations, we consider each system to be excited by an 
incident plane wave, as represented in the following scheme. 

Oo=0, a£(n, L) = 8 })p >8 nn >8(L — L') ; 

ao = 0, a p {n, L) =8 pp »8 n n»8(l> — L ;/ ). 

From (9) we find 

b p ,,(n'\ — h n )y) H »(L") =b"j ) ,{n\ — \u')j]n'{L') 

and from (10) and (11) 

b' q (m, K)=Sq P '(m, K; ai', L'), b" q {m, K) =S qp >'(m, K; //", L"). 

Hence (after changing variables to get rid of the primes), we obtain the set of scattering reciprocity 
relations 

i ?m («)S pg (m,K;ii,L) = 7j JI (L)S&(/i,-L;nt,-K). (A-13) 

In this case interchanging S and S a yields no further information. 

If the reader has followed through any one of the above exercises, he will have no trouble in 
showing 

Soo = Sft. (A-14) 

If we introduce the reciprocity dyadic 

tl = 7}iK\Ki + T)iK2K2 (A-15) 
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and use the definitions (1.3-3), then (for example) (12b) becomes 

Tj So,(K) = i|-S«o(-K) (A-16) 

and (13) is summed up in 

rf(W) •S, / (k,L) = S;(-L,-K)^r 1 (L), (A-17) 

where ij(K) denotes 17 as a function of K and the superscript 'T" denotes the transposed dyadic 
(obtained by transposing the elementary dyads involved). It is worth noting that the minus sign 
appearing between the members of (12b) does not appear in (16). The reciprocity dyadic is 
manifestly diagonal in the K\ , k> basis; (4.2-3) shows what happens when r\ is presented in the 
e x , Gy basis. 

The Lorentz reciprocity relation seems to have been the quite generally preferred basis for 
the derivation of transmission or coupling equations (similar, at least in function, to our (2-15)) 
[19, 21, 25, 26]. There is, however, no good reason why reciprocity should be invoked for that 
purpose: receiving characteristics can be defined analytically and operationally, independent of 
transmitting characteristics (as is done in this paper). Nevertheless, it is of some interest that the 
generalized Lorentz relation (8), or, more conveniently in the plane-wave framework, the lemma (9), 
can be used to derive transmission equations without invoking (ordinary) reciprocity. 

5.2. Appendix B: Dissipative Characteristics of Media; Comparison for Mutually Adjoint Media 

The concept that a medium may be "lossy," "lossless," or even "gainy" at a given point is 
familiar. We first need to give the mathematical expression of these qualitative properties of 
media. Then it will be easy to show that these properties are point-wise identical for mutually 
adjoint media. 

According to properties of the complex Poynting's vector, dissipation or power loss per unit 
volume, #(r), is given by 

g(r)=-|Re V • (ExH). (B-l) 

In the medium described by (1.5-1), Maxwell's equations require 

V x E =ico(ti • H + i> • E), 

(B-2) 

V x H = -io)(€ E + r H). 

(Note that these equations are source-free in the sense that they are satisfied by E = H = 0.) 
Using a standard vector identity, we combine (1) and (2) to obtain 

q(r) = -^Re {ico[H -|t-H-E-c-E + H- (*>-t*) • E]}, (B-3) 

where the superscript "*" denotes the Hermitian conjugate. For our purpose it is convenient to 
rearrange this expression to read 
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The dissipative properties of the medium are thus seen to be determined by properties of the 
Hermitian matrix 



A(r) 



V — 






V * — T 



(B-5) 



which we call the loss matrix for the medium. For the present purpose the characteristics of the 
medium at a given point are appropriately classified according to the values assumable by q(r) 
as E and H, considered independent and arbitrary, vary. All possibilities are listed and named 
in table B-l. 

There is no well-established terminology precisely fitting the physical properties being dis- 
cussed. Terms generally chosen, and indeed the term "loss matrix," are biased by the tacit assump- 
tion that a medium normally is lossy. Thus if the loss matrix happens to be the zero matrix, we 
would ordinarily say that the medium is lossless at the point in question. But the medium would 
also be gainless at the same point. Neither physical nor technological restrictions prevent the 
occurrence of negative energy loss, which might be called gain (cf. Brand [51]). 

Mathematical criteria for deciding the value class of Hermitian matrices are given by Mirsky 
[52]. In the lossless case, which is physically important but not counted as defining a value class 
by Mirsky, the loss matrix must be the zero matrix; this in turn requires € = €*, t= i>*, and /£ = /**. 



Table B-l. Classification of dissipative characteristics 



Values assumed 


Value class of loss 


Dissipative characteristics 


by </(r)„ 
r fixed 




matrix at point r 


of medium at point r 


>0 




Positive definite 


Unconditionally | 


^0 




Positive semidefinite 


> lossy or passive 
Conditionally J 


= 




Zero 


Lossless, gainless or neutral 


^0 




Negative semidefinite 


Conditionally 


<0 




Negative definite 


j- gainy or active 
Unconditionally ) 


>0 
< 




Indefinite 


Indefinite 



Note: In the semidefinite cases both the zero and the nonzero values are to be as- 
sumable, and in the indefinite case both positive and negative values are to be assumable. 

The medium adjoint to that described by (1.5-1) is defined and described by (1.5-2), according 
to which the adjoint expressions can be obtained by the replacements €~ >€, t^~v, i>— » — t, 
and /*- > jt*. For the loss matrix this yields 



A«(r) = A(r) 



(B-6) 



i.e., A"(r) is the transpose of A(r). This relation is sufficient to insure that the corresponding loss 
functions, q(r) and q a (r), belong to the same value class (this can be seen, e.g., from the criteria 
given in Mirsky). We may say that the dissipative characteristics of a medium and its adjoint are 
the same, point by point. Whatever distribution of characteristics, including regions of active 
media (as in some antennas), is realized in one system will be realized in the adjoint system. 

5.3. Appendix C: Two-Dimensional, Spatial Sampling or Interpolation Theorem 

A variety of two-dimensional sampling theorems may be found in the literature [42, 53]. The 
methods of derivation often used are unnecessarily complicated for our purposes. For the 
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convenience of the reader we sketch a simple derivation of the simple desired result. 

With the abbreviation S ( ' )2 (K) • Sio(K) exp (iyd) =/(K), the equation under discussion. 
(3.1-1), becomes 



a F' ff(K 



6i(P) = a F ' j /(K) e' K 'V/K. (C-l ) 

We have already noted that the mathematical requirement that /(K)ibe "band limited" may be 
fulfilled extraordinarily and almost arbitrarily well. Let us therefore assume that nonzero values 
of /(K) occur only in a finite region K of wavenumber space. For simplicity, we take K to be 
rectangular, bounded by the lines /.., — ±k\ , k y = ±/r 2 . Then (1) may be written 

b' Q (P) = a F ' j' 1 j kl /(K) e' K • *dk x dk y . (C-2) 

Further, /(K) presumptively can be represented in K by a double Fourier series, with periods 
2/. i and 2/»2 : 

/(K)= J J Crse ,k ,Vs - ,(: ~ 3) 

r = — x x = — x 

where the coefficients are given by the usual formula 

C rs 



I/,,/. 
/ith 



- f' 2 r /(K)c« k ' p « dk x dky (C-4) 

2 J A- a J-A, 



P„ = p', + pv (C-5) 

Comparing (4) and (2), we see that the C r . s are proportional to 6 '(P) evaluated for P=P rs , so 
that (3) may be written 

/(K)= / ,., £ b o(Pr S )e -«'K-p„ (C-6) 

Thus/(K) is completely determined by the sampling of data at the discrete points P rs - This is 
the result used in (3.1-4). 

The usual objective in sampling theory is the reconstruction of the sampled function. In the 
present instance this is easily accomplished by substituting (6) in (2) and integrating. The result 
is an expression involving cardinal functions and expressing b' (P) at all points in terms of its values 
at the lattice points. 

One should observe the reciprocal relation between the size of the region K in wavenumber 
space and the size of the elementary cell in .r, y space. The smaller the former, the larger the lat- 
ter—and the smaller the required density of sampling points in x, y space. 
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